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OBJECTIVE 


Assess  the  feasibility,  and  then  the  utility,  of  inverting  the  far  field  patterns  of 
shipboard  antennas. 

DISCUSSION 

Prony  processing  may  be  applied  to  antenna  patterns  to  determine  source  locations 
and  can  accurately  identify  the  principal  radiators  even  for  the  cluttered  shipboard  environ- 
ment. The  technique  could  be  a very  useful  tool  during  the  shipboard  antenna  arrangement 
design  in  cases  where  one  or  more  parasitic  structures  are  Rresent.  Not  only  can  these 
structures  be  located  by  the  Prony  poles,  but  their  relative  importance  can  be  deduced 
from  the  size  of  the  residues. 

Prony  processing  is  not  the  only  possible  approach,  as  demonstrated  via  Fourier 
series  processing.  But  the  chief  advantages  of  Prony  are  that  (1 ) unlike  the  Fourier  series 
approach,  Prony  processing  can  resolve  sources  separated  by  less  than  one-half  wavelength, 
and  (2)  the  solution  is  found  with  fewer  iterations.  Perhaps  with  experience  or  some 
suitable  eigenvalue  technique,  as  suggested  by  Lawrence  Livermore  Laboratory,  only  one 
or  two  iterations  may  be  required.  The  only  disadvantage  of  Prony  is  its  susceptibility  to 
noise  in  the  pattern  data.  This  problem  is  now  receiving  considerable  attention  from  the 
technical  community. 
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1.0  INTRODUCTION 


1.1  OVERVIEW 

The  far  field  pattern  of  an  electromagnetic  source  is  a Fourier  transform  of  the 
spatial  distribution  of  the  source.  * Because  patterns  can  be  expressed  in  terms  of  an 
exponential  series  of  the  spatial  distribution,  they  are  candidates  for  data  transformation 
or  “inversion.”  An  inversion  process  which  has  been  found  to  be  increasingly  useful  is 
one  based  upon  Prony’s  method.  This  procedure  is  one  by  which  the  coefficients  and 
exponents  of  an  exponential  series  can  be  derived  from  a sequence  of  equally  spaced 
sampled  values. 

Data  inversion  could  be  useful  in  two  ways:  data  characterization  a>’i  data  com- 
pression.” Prony  processing  can  be  used  to  characterize  the  source  distribution  of  a far 
field  pattern.  Dr.  E.  K.  Miller  has  indicated  that  the  exponents  in  the  spatial  domain 
correspond  to  the  physical  locations  of  the  sources  in  real  spaced  One  application  might 
be  the  imaging  of  a radiation  source.  Another  possibility  is  that  of  developing  a model 
for  the  object  that  in  a perturbation  sense  could  be  used  to  predict  its  new  behavior  when 
its  geometry  is  modified.  A specific  application  for  this  capability  would  be  the  use  of 
data  obtained  on  the  NOSC  scale  model  range,  to  design  and  evaluate  antenna  systems 
on  ships.^ 

The  Prony’s  process  might  also  be  used  to  compress  the  original  data.  In  this  case, 
the  coefficients  and  exponents  may  have  no  particular  physical  meaning,  but  would  repre- 
sent a minimum  amount  of  numbers  required  to  store  and  manipulate  the  original  data. 
Thus  ship’s  pattern  data  could  be  stored  in  a considerably  reduced  quantity,  eg,  instead 
of  storing  360  complex  numbers  per  pattern,  it  may  be  necessary  to  store  only  4 to  10 
pairs  of  complex  Prony  coefficients.  This  application  could  result  in  more  efficient 
antenna  optimization  or  direction  finding  algorithms. 

1.2  OBJECTIVE 

The  technical  objective  of  this  study  is  to  assess  the  feasibility  and  then  the  utility 
of  inverting  the  far  field  patterns  of  shipboard  antennas.  Two  applications  are  considered: 
data  characterization  and  data  compression. 


1 . Kraus.  JD,  Antennas.  McGraw  Hill  Book  Company,  Inc..  1950. 

2.  Lawrence  Livermore  Laboratory  Report,  UCID-1751 1,  Data  Characterization  and  Compression, 
by  EK  Miller.  July  6.  1977. 

3.  Lawrence  Livermore  Laboratory  Report,  UCID-1750s,  Prony’s  Method  for  Angle  Domain,  by 
tK  Miller  and  DL  Lager,  June  14,  1977. 

4.  Rockway.  JW  and  DuBrul,  DW.  Performance  Prediction  Analysis  for  Shipboard  Antenna  Systems. 
Naval  Engineers  Journal,  October  1977. 
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1.3  APPROACH 

A number  of  alternatives  to  Prony  processing  of  pattern  data  exist.  However, 
evaluation  and  comparison  of  Prony-type  processes  is  beyond  the  scope  of  this  study, 
except  where  a similar  or  related  technique  might  provide  insight  into  the  Prony  process. 

An  example  is  the  evaluation  of  pattern  data  using  a truncated  Fourier  series.  This 
technique  is  considered  along  with  the  Prony  method.  Dr.  Miller  has  suggested  that  the 
Fourier  series  approach  is  a special  case  of  the  more  general  Prony  process.^  Thus,  the 
more  familiar  and  extensive  literature  on  the  Fourier  series  can  be  used  to  provide  in- 
sight into  the  Prony  process.  Both  the  Fourier  series  and  the  Prony  methods  are  discussed 
in  detail  in  section  2.  Both  methods  have  been  implemented  into  a computer  code  called 
Angle  Domain  Pattern  Analysis  Program  (ADPAP).  This  code  is  also  discussed  in  section  2. 

A user’s  guide  and  code  listing  for  ADPAP  are  given  in  appendices  B and  C respectively. 

Two  applications  of  antenna  pattern  data  inversion  are  considered:  data  character- 
ization and  data  compression.  Section  3 presents  the  application  to  pattern  data  character- 
ization. Section  4 presents  the  application  to  pattern  data  compression.  In  section  5,  the 
feasibility  and  then  the  utility  of  pattern  data  inversion  are  summarized. 

2.0  DATA  INVERSION  METHODS 

2.1  PRONY’S  METHOD 

A ship  is  an  electrically  long  narrow  conducting  platform  with  many  definite  antenna 
and  parasitic  protrusiles.  A single  transmitter  driven  antenna  may  excite  one  or  more  para- 
sites. The  result  is  an  array  which  is  nearly  collinear.  As  a first  approximation,  the  shipboard 
antenna  problem  will  be  treated  using  elementary  array  theory. 

The  far  field  of  a collinear  array  of  n isotropic  point  sources  having  coordinates  XQ 
and  complex  amplitudes  AQ.  for  a = 0,  1,2...,  n-1  is 

E(0)  = | A e'0cos0 
a=0 

where  d = 2ir/\  and  0 is  the  observation  angle  relative  to  the  array  axis.  The  pattern  inversion 
problem  is  to  determine  the  locations.  Xa.  and  the  eomples  amplitudes,  AQ1.  of  an  array  of 
sources  which  produce  the  far  field  L(0). 

Prony ’s  method  is  a technique  for  approximating  a set  of  data  as  an  exponential 
series.  First  the  data  must  exist  in  equally  spaced  samples  of  cos  0.  A linear  interpolation 
scheme  is  sufficient  to  convert  from  equal  spacing  in  0 to  equal  spacing  in  cos  0.  Taking 
advantage  of  the  symmetrical  nature  of  the  patterns  to  avoid  multiple  values  of  cos  6,  the 
data  is  restricted  to  the  interval  0 < 0 < n. 

Let  6 be  the  increment  of  cos  0 for  N samples  and,  for  convenience,  let  qj.  = e*^xk^ 
where  k = a + 1 . Then  ( I ) becomes  a system  of  equations: 


5.  Lawrence  Livermore  Laboratory . Private  Communication  by  EK  Miller.  June  1978. 
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Eo  = 2 Ak 

k=l 


E]  = I Akqk 

k=l 


E2  = I Akqk 

k=l 


En-1  _ 2 Akqk  1 
k=l 

If  the  q’s  were  known,  (2)  would  be  a system  of  linear  equations  and  the  solution  straight- 
forward. However,  in  this  case,  the  q’s  are  unknown  and  (2)  is  nonlinear  in  q. 

Next,  assume  that  q j , 1 2 ■ . ■ qn  are  the  roots  of  the  equation: 

qn  - a lqn~l  - a2qn-2  - . . . -an_,q  - an  = 0 (3) 

Now  multiply  the  first  equation  of  (2)  by  -o:n,  the  second  by  -an_j. . . ■ the  ( n—  1 )th  by 
-ai,  and  the  nth  by  +1  and  add  the  results.  Using  the  relationship  of  (3).  the  result  is: 

En  = alEn-l  + + anEo  = akVk  (4) 

k=l 

Similarly,  proceed  starting  with  the  second,  third,  fourth,  etc.,  to  the  (N-n)  th  equation. 
The  result  is  a linear  system  of  equations: 


En  2 akEn-k 
k=l 


En+1_  2 akEn+l-k 

k=l 


eN-1  ".I  akEN-l-k 


Equations  (5)  are  linear  in  a and  are  easily  solved.  The  values  of  q are  then  readily  deter- 
mined as  the  roots  of  (3).  The  source  locations  Xa,  also  called  poles,  are  found  from 

qk  = e^k5 

where  k = a + 1 . Equation  (2)  can  then  be  solved  for  the  complex  amplitudes  Ak,  also  called 
residues,  of  ( 1 ).  Reference  6 gives  a more  complete  discussion  of  Prony’s  method.  Both  the 
poles  and  residues  are  complex.  The  actual  sources  will  have  higher  residues  and  purely 
imaginary  poles.  The  other  poles  and  residues  which  may  appear  are  merely  curve  fitting 
values  and  have  no  physical  interpretation. 

2.2  FOURIER  SERIES  METHOD 

Analysis  of  pattern  data  using  a truncated  Fourier  series  is  an  alternative  approach  to 
identifying  the  source  locations.7  The  approach  is  somewhat  simpler  than  Prony  ^nd  provides 
different  insight  into  pattern  inversion. 

The  far  field  of  a collinear  array  of  n isotropic  sources  is  given  in  ( 1 ).  Let  the  ampli- 
tude be  given  by  Aa  = a^e'0"  where  </>a  is  the  phase  of  the  source  a.  Then  ( 1 ) becomes: 

n">  i(/3Xacos  e + <t>a) 

E(0)=  2 aQe  a “ (6) 

o=  0 

Noting  that  E(0)  = Er(0)  + i E|  (0),  the  real  and  imaginary  parts  of  the  field  can  be  expanded 
separately: 


Er(0)  = £ aa  cos  (|3Xa  cos  0 + 0a) 
ck=0 


£ [(aa  cos  0Q)  cos  (0Xa  cos  0) 
a=0 

Haa  sin  da)  sin  (/3Xa  cos  0)1 


Letting 


•o  = aocos0o 


Ja  = -a  sin  0a 


6.  Hildebrand.  FB,  Introduction  to  Numerical  Analysis.  McGraw-Hill  Book  Company,  Inc.,  1956. 

7.  Kreyszig,  E.,  Advanced  Engineering  Mathematics.  John  Wiley  and  Sons.  Inc.,  1962. 


(7)  becomes 


n-i 

Er(0)  = ^ [Ia  cos  (/3Xa  cos  0)  + JQ  sin  (/3Xa  cos  0)]  (9) 

a=0 

Assuming  the  sources  are  equally  spaced,  let  Xa  = ad  and  u = /3d  cos  0,  then  (9)  becomes: 

n-1  n-1 

Er(0)  = IQ+  X Iacos(au)+  £ JQ  sin  (au)  (10) 

a= 1 a= 1 

Equation  ( 10)  is  a truncated  Fourier  series  where  the  coefficients  can  be  found  in  the  usual 
way: 


I 


o 


7T 

J Er(u) du 
—it 


1 


a 


/E  r<“> 


cos (au) du 


-7T 


J 


a 


u)  sin  (au)  du 


(11) 


The  spacing  of  the  sources  and  the  angular  spread  of  the  pattern  are  critical  to  the  evalu- 
ation of  ( 1 1 ).  Taking  advantage  of  the  symmetrical  nature  of  the  patterns,  the  pattern 
data  may  be  restricted  to  the  interval  0 < 0 < 7r.  The  minimum  source  spacing  for  (11) 
to  be  valid  is  d = A/2.  This  is  the  minimum  separation  between  sources  that  can  be 
physically  resolved  based  on  patterns. 

The  integrals  in  (1 1)  can  be  evaluated  using  numerical  integration.  First  subdivide 
the  interval  of  integration  into  q equal  parts  of  length  6u  = 7t(6  cos  0).  The  sampling  of 
the  pattern  data  must  be  in  equal  increments  of  cos  0 as  in  the  Prony  process.  The  sample 

points  are  located  at  uQ  = n cos  (o° ),  u j , u2 uq  = 7r  cos  (rr).  The  first  integral  in  ( 1 1 ) 

is  found  by: 

( 6 cos  0 ) 

lG  = lEr(u0)  + Er(uj)  + ....Er(uq)] 

similarly 

la  = (6  cos 0)  [Er (uQ) cos au0 + . . . . Er (uq) cos auq] 

Ja  = (5  cos  0)  [Er  (uQ)  sin  auQ  + 


Er  (uq)  sin  auq) 


Associated  with  each  pair  of  coefficients,  (la,  Ja)  is  a source  location. 


xa  = ±a(X/2) 

(12) 

The  source  amplitude  and  phase  are 

aa  =<IS  + Ja)1/2 

(13) 

<pa  = tan-1 

( 14) 

A similar  set  of  equations  can  be  developed  for  the  imaginary  part  of  the  pattern. 
n-1 

E[  (0)  = X aa  sin  (^Xa  cos  6 + 
a=0 

n-1 

= £ [(aa  sin  </>a)  cos  (0Xa  cos  0)  (15) 

a=0 

+ (aa  cos  <pa)  sin  (PXa  cos0)| 

Assume 

la  = aa  sin  <pa 


•*a  aa  cos 
Xa  =ad 

U = 7T  COS  0 


with  d = — and  0 < 0 < 7T. 

L 

The  result  is: 


n-1  n-1 

E|(u)  =l0+£  cos  (au)  + £ Ja  sin(au) 
a=  1 of  1 


(16) 


where  the  coefficients,  0a.  JQ)  are  evaluated  as  before.  Again,  a source  location  is 
associated  with  each  pair  (la.  Ja)  with  amplitude  and  phase  of 
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a = n2  + J2>V4 
aa  ‘‘a  Jcr 

4>a  = tan”1  (la/ia)  (17) 


An  ambiguity  in  the  source  locations  exists  if  this  approach  is  directly  applied  to 
pattern  inversion.  Equation  ( 1 2)  indicates  the  problem.  Two  source  locations  are  actually 
associated  with  each  pair  of  coefficients.  What  may  occur  is  that  two  different  arrays 
having  different  source  locations  and  source  strengths  can  produce  identical  patterns.  The 
difficulty  is  easily  averted  by  premultiplying  the  pattern  data  by  a phase  shift  e'^s,  which 
in  effect  shifts  the  pattern  phase  center  far  enough  along  the  assumed  array  axis  so  that 
the  true  source  locations  are  revealed.  Of  course,  this  requires  some  prior  knowledge  of 
the  spatial  extent  of  the  actual  source  so  that  the  shift  is  sufficient. 

2.3  ANGLE  DOMAIN  PATTERN  ANALYSIS  PROGRAM  (ADPAP) 

The  Prony  method  of  section  2.1  and  the  Fourier  Series  method  of  section  2.2  have 
been  coded  in  the  Angle  Domain  Pattern  Analysis  Program  (ADPAP).  A user’s  guide  and  pro- 
gram listing  for  ADPAP  are  given  in  appendices  A and  B,  respectively.  The  Prony  analysis 
portion  of  ADPAP  is  adapted  from  the  Lawrence  Livermore  Laboratory  code  SEMPEX.^ 

The  Fourier  series  capability  and  shifting  theorem  coding  were  written  at  NOSC  using  the 
approach  of  section  2.2. 

3.0  APPLICATION  TO  PATTERN  DATA  CHARACTERIZATION 

The  NOSC  scale  model  pattern  range  has  the  capability  to  measure  both  magnitude 
and  phase  of  the  scale  model  antenna  patterns.  However,  to  date,  pattern  data  with  both 
magnitude  and  phase  information  does  not  exist.  Also,  time  prohibited  collection  of  such 
data  for  a wide  range  of  antenna  configurations. 

Q 

Fortunately,  a method  of  moments  computer  code,  NEC  . is  available  on  the  NOSC 
computer.  NEC  has  equivalent  or  better  accuracy  than  the  scale  model  range  and  the  ver- 
satility to  model  a wide  range  of  antenna  configurations  in  a timely  manner.  Consequently, 
NEC  was  used  as  the  source  of  pattern  data  for  the  evaluation  of  the  Prony  and  Fourier 
techniques. 

3.1  SIMPLIFIED  ARRAYS 

The  first  problems  considered  were  simple  two-  and  three-element  arrays.  As 
confidence  and  insight  were  developed,  more  complex  configurations  were  considered.  A 
few  of  these  elementary  arrays  (figures  1 , 2,  and  3)  have  been  selected  for  illustration. 


8.  Air  Force  Weapons  Laboratory,  Kirtland  Air  Force  Base,  New  Mexico.  AFWL-TR-76-200.  Users  Manual 
for  SEMPEX:  A Computer  Code  for  Extracting  Complex  Exponentials  from  a Time  Wave  Form,  by 
DL  Lager,  HG  Hudson,  and  AJ  Poggio,  1976 

9.  Naval  Ocean  Systems  Center  Technical  Document  1 16,  Numerical  Electromagnetic  Code  (NEC)  Method 
of  Moments,  by  GJ  Burke  and  AJ  Poggio.  18  July  1977. 
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Figure  1 shows  two  half-wave  dipoles  separated  by  three  wavelengths  in  free  space. 

The  dipoles  are  center-fed  in  phase.  The  far  field  pattern  in  the  X-Y  plane,  the  magnitude 
and  phase  of  the  vertical  component  of  the  electric  field,  was  calculated  using  NEC.  Then 
the  Prony  process  and  Fourier  series  analysis  were  applied  to  the  pattern  data  using  ADPAP. 

The  results  are  shown  in  tables  1 and  2. 

• The  residues  and  poles  listed  in  table  1 are  related  to  the  source  patterns  and  locations 

by  equation  ( 1 ).  Residues  are  the  amplitudes,  Aa,  and  the  poles  are  the  arguments  of  the 
exponentials.  The  residues  correspond  to  the  magnitude  and  phase  of  possible  sources,  and 
the  poles  correspond  to  the  source  locations  in  meters  in  a complex  spatial  domain.  Real 
sources  have  poles  on  the  real  axis.  In  this  case,  pole  numbers  1 and  2 have  the  highest 
relative  residues  and  are  therefore  prime  source  candidates.  The  real  parts  of  these  poles  are 
nearly  zero  compared  to  the  other  poles  in  the  list  The  imaginary  part  gives  the  source 
location  as  ± 1.5  meters  along  the  array  axis.  Note  that  the  source  strengths  are  the  same 
magnitude  and  phase.  The  other  poles  listed  appear  as  curve  fitting  poles  and  have  no  re- 
lationship to  physical  sources. 

The  solution  of  the  matrices  in  the  Prony  process  used  by  ADPAP  is  presently  set 
up  for  square  matrices.  This  requires  that  the  number  of  pattern  data  samples  be  twice  the 
number  of  poles  specified.  Errors  are  introduced  by  specifying  too  many  or  too  few  poles, 
either  because  of  insufficient  data  samples  to  characterize  the  pattern  or  because  over- 
specification of  poles  introduces  curve  fitting  noise.  Without  some  a priori  knowledge  of 
the  number  of  poles,  and  hence  data  samples  required,  the  Prony  process  must  be  applied 
repeatedly  until  a stable  solution  can  be  obtained.  Of  course  this  requires  some  knowledge 
of  the  source  distribution  in  the  first  place,  but  this  is  an  existing  condition  for  this  par- 
ticular application.  Table  1 is  the  result  of  several  iterations  of  the  Prony  process  applied 
to  the  pattern  data  by  specifying  a different  number  of  poles  (eg,  number  of  poles  = 10, 

15,  20,  25, . . .). 

The  Fourier  coefficients  listed  in  table  2 are  related  to  the  pattern  data  of  the  array 
in  figure  1 by  equation  (10).  They  were  found  by  expanding  the  real  component  of  the 
electric  field  in  a Fourier  series  as  discussed  in  section  2.2.  The  source  strength  and  location 
associated  with  each  pair  of  Fourier  coefficients,  1(N)  and  J(N),  are  given  by  equations  (12), 
(13),  and  (14).  The  values  of  location  and  strength  are  listed  in  the  columns  under  S,  AMP  (N), 
and  PHASE,  respectively.  A different  set  of  Fourier  coefficients  could  also  be  obtained  by 
expanding  the  imaginary  component  but  the  source  strengths  and  locations  are  identical. 

Hence  the  examples  discussed  will  use  the  Fourier  series  expansion  of  the  real  component 
only. 

One  of  the  options  of  ADPAP  is  to  apply  a shift  of  exp  (i/3X  cos  0 ) to  the  pattern 
» in  n equal  steps  prior  to  computation  of  the  Fourier  coefficients.  In  effect,  if  the  observer 

is  on  the  array  axis  (the  X axis  in  figure  1 ),  the  shift  moves  the  phase  center  away  from  the 
observer.  The  source  locations,  S,  in  table  2 are  given  with  respect  to  the  original  phase 
center.  Five  sets  of  coefficients  are  shown  for  varying  amounts  of  shift  from  zero  to  2.0 
meters.  Originally  the  shift  was  applied  in  equal  increments  of  0.1  meters,  but  in  the 
, 1 interest  of  conserving  space,  only  five  of  the  steps  are  illustrated. 

The  Fourier  coefficients,  and  the  source  amplitudes  as  a result,  rise  and  fall  as  the 
phase  center  is  shifted  along  the  array  axis.  Larger  amplitudes  indicate  possible  source 
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Table  1.  Residues  and  poles  from  Prony  processing  of  the  pattern  of  two  half-wave  dipoles. 
(See  figure  1 .) 


COMPLEX  PRONY 


FRECUENCY  = 

.2998+009 

HZ 

WAVELENGTH  = 

1 .0000 

*et  e*  s 

NO . OF  SAMPLE S = 

40 

NO.  Of  POLES  = 

2 G 

RESIDUES 

°0L  E S 

MAG 

PHASE 

real 

I M A G 

1 

.59669+002 

89.794 

.8  6 44  7 -003 

.1 '006+001 

2 

.59676+002 

90.205 

.“6118-003 

-.1 '006+001 

3 

.1  29  6 2*001 

94.039 

-.10942+001 

-.62827+001 

4 

.1 C098+001 

-41.821 

-.1  369‘7  + 001 

.61389+001 

5 

.23319+001 

129.81  3 

-.1  0832  + 001 

-.19335+001 

6 

.2  tB62  + 0Q1 

-43.000 

-.1  0993+001 

.54456+000 

7 

.14580+001 

-157.482 

-.1322«+001 

.41 318+001 

8 

.10373+001 

-83.162 

-.14082+001 

-.51536+001 

9 

.1  1945  + 001 

14  7.886 

-.97014+000 

.93118+001 

10 

.1  1928+001 

32.136 

-.969S«  +C00 

-.931 16+001 

1 1 

.28861+001 

-137.003 

-.10991+001 

-.5451 8*000 

1 2 

.23323+001 

50.163 

-.1  0?9"+0C1 

.19325+001 

1 3 

.1 9086+001 

56.176 

-.12*63+001 

-.39200+001 

14 

.1 91Q5+001 

123.833 

-.1  2360  + 001 

.3ni 94+001 

1 5 

.14580+001 

-22.495 

-.13228+001 

-.41 324+001 

16 

.10383+001 

-96.833 

-.14081+001 

.51525+001 

17 

.1 0084+001 

-135.131 

-.13700+001 

-.61395+001 

18 

.1  1563  + 001 

f 8.633 

-.12416+001 

.71772+001 

19 

.1  1551  + 001 

161.450 

- . 1 2 4 1 7 ♦ 001 

-.71775+001 

20 

.1 2974+0C1 

86.05 1 

-.10944+001 

.87“ 25+001 
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Table  2.  Fourier  coefficients  for  the  pattern  of  two  half-wave  dipoles,  (see  figure  1 .) 


FOURIER  COEFFICIENTS  FOP  R E AL  ( E ) 

FREQUENCY  - .2998+009  HZ 

WAVELENGTH  = 1.0000  METERS 

NO.  OF  SAMPLES  = 100 

NO.  OF  TERMS  = 10 


The  phase  center  HAS  BEEN  SHIFTED  by  .OCGO  meters  along  array  axis 


N 

S(METERS) 

1 (N) 

J (N) 

A*P(N) 

PHASE 

0 

.0000 

-.51945-004 

.1  21  59  + 001 

-90.002 

1 

.5000 

.16184-005 

.30897-011 

.1 61 84-005 

-.000 

2 

1 .0000 

-.16706-005 

.54534-012 

.1 6706-005 

-180.000 

3 

1 .5000 

. 1 6205-005 

.1  1 826-011 

.1  6205-005 

.000 

4 

2 .0000 

-.15867-005 

- .99932-012 

.1 5R67-C05 

180  ."DO 

5 

2. <000 

.16922-005 

.63953-011 

.1  6922-005 

-.000 

6 

3.0000 

-.3555 =-004 

.61  097-011 

.3555®- 004 

-180.000 

7 

3. 5000 

.16931-005 

.1  4 584-01  1 

.1  6931-005 

.000 

3 

4.0000 

-.15937-005 

-.55  356-01  1 

.1  5937-005 

180.000 

9 

4.5000 

.16519-005 

.22368-011 

.16519-005 

.000 

10 

5 .0000 

-.16653-005 

- .4862  5-011 

.1  6653-005 

180.000 

* Hit  * it  X-  ********* 


THE 

PHASE  CENTER 

HAS  BEEN  SHIFTED  BY  .5000 

METERS  along 

ARRAY  AXIS 

N 

S (METE  RS  ) 

I (N) 

J < \) 

AhPCn  ) 

PHASE 

0 

1 

-:o8S8 

:;3USi:88I 

-.2093  1-001 

.1 2936+001 
.20931-001 

90.000 

90.140 

2 

.5000 

.9  3947-005 

.59906+002 

.62630+002 

-90.000 

3 

4 

!:?888 

.35  882-006 
-.75523-005 

-.71930-002 
- .59904+002 

89.997 

90.000 

5 

2.0000 

-.20108-004 

.2856  8-001 

.2 856«-0G1 

-90.040 

6 

2.5000 

.14019-005 

- .53875-002 

.53881-002 

89.984 

7 

3.0000 

-.16921-004 

.40430-003 

.4  0465-003 

-92.397 

8 

3. 5000 

.12258-005 

- .1  3001-002 

.1 3001-002 

89.946 

9 

4.0000 

-.10633-005 

.21 349-002 

.21348-002 

-90.029 

10 

4.5000 

.71527-006 

- .64536-002 

.64534-002 

89.994 

*********** 

* * • • 

r 

- — — - in  ..jjjpgg 

Table  2 

Continued. 

j 

THE 

PHASE  CENTER 

HAS  BEEN  SHIFTED 

BY  1.5000 

KETERS  AL0N6 

ARRAY  axis 

N 

S(METERS) 

I (N) 

J (N) 

a*p(n) 

PHASE 

0 

-1 .5000 

.75740-005 

.63011*002 

90.000 

1 

-1.0000 

.25672-005 

-.71  73  7-002 

.71736-002 

89.979 

2 

-.5000 

-.14016-005 

.1  7047-002 

.1 7047-002 

-90.047 

3 

.0000 

-.68424-004 

.4435  7-003 

.44882-003 

-98.769 

4 

.5000 

.11080-006 

- .36705-002 

.36701-002 

89  .99» 

5 

1.0000 

.39012-006 

.21  76  5-001 

.21761-001 

-89.999 

6 

1.5000 

-.80195-005 

- .5991  1 ♦ 002 

.60243*002 

90.000 

7 

2.0000 

-.13342-005 

.31 11 9-001 

.31122-001 

-90 .002 

8 

2.5000 

-.39975-006 

- .1  3150-001 

.1 31 51-C01 

90 .002 

9 

3.0000 

-.12080-004 

- .22106-002 

.221 06-002 

90.313 

10 

3.5000 

-.53943-006 

-.1  1006-001 

.1  1007-001 

90 .00'* 

* ************** 

THE 

PHASE  CENTER 

HAS  BEEN  SHIFTED 

BY  1.6000 

WETERS  along 

array  axis 

N 

S(WETERS) 

I (N) 

J ( N) 

AHP ( N ) 

PHASE 

0 

-1  .6000 

.38147-005 

.62952*002 

90.000 

1 

-1.1000 

.93055-005 

- .2274  8*002 

23072*002 

90.000 

2 

-.6000 

-.66661-005 

.10004*002 

.99559*001 

-90.000 

3 

-.1000 

-.54735-004 

- .51  630*001 

.51813*001 

90.001 

4 

.4000 

-.18804-004 

.1561 8*001 

.1  5619*001 

-90.001 

5 

.9000 

.14527-004 

.39396*001 

. 3 94  58 ♦ 001 

-90.000 

6 

1 .4000 

-.18977-004 

- .51470*002 

.51077*002 

90.000 

7 

1.9000 

.63193-005 

- .1  7939*002 

.1  7621  * 002 

90 .0Cn 

8 

2.4000 

-.78799-005 

.96998*301 

.96263*001 

-90.000 

9 

2.9000 

-.26973-006 

- .71035*001 

.61706*  001 

90.000 

10 

3.4000 

-.15875-004 

.57231*001 

.57299*001 

-90.000 

the  phase  center  has  been  shifted  by  2.0000 

* * 

K E T E & S A L ON  6 ARRAY  AXIS 

N 

S(METERS) 

I (N) 

J ( N) 

AWP(N  ) 

phase 

0 

-2.0000 

.15991-005 

.1  2469*  001 

-90.000 

1 

-1.5000 

.65409-005 

- .59904*00 2 

.62117*002 

90.000 

2 

-1.0000 

-.15503-004 

.21  39  0-001 

.21390-001 

-90.04? 

3 

-.5000 

-.97475-006 

-.36309-002 

.36®  08-002 

90.014 

4 

.0000 

-.49394-004 

.84  596-003 

.34741-003 

-9  3 .34  ? 

5 

.5000 

-.6261 2-006 

- .4969  3-007 

.49691-002 

90.007 

6 

1.0000 

-.6261 ® -00  6 

.2  T 6 9 7-00 1 

.23391  -001 

-90.00? 

7 

1 .5000 

-.10274-004 

- .5991  7*002 

.5  7 7 75  ♦ C 02 

90.000 

8 

2.0000 

. T9399-005 

.26372-001 

.26371  -001 

-69.991 

9 

2.5000 

-.3515® -006 

-.1  6 399-001 

.1  6397-001 

90.001 

10 

3.0000 

-.16121-004 

- .1  1 239-001 

.1  1240-001 

90.08? 

^ — 

..  ..  : 

12 

r — 

locations.  By  stepping  in  fractional  wavelengths  through  a total  shift  of  at  least  one-half 
wavelength,  the  source  locations  can  be  determined  within  an  error  equal  to  the  step  size. 

A maximum  in  amplitude  occurs  when  the  coefficients  coincide  with  real  source  locations. 

The  first  set  of  coefficients  in  table  2 are  the  results  for  zero  shift.  The  first  co- 
efficient has  a significant  amplitude  compared  to  the  amplitudes  of  the  other  coefficients 
indicating  that  a single  source  could  be  located  at  the  origin.  However,  as  a shift  is  applied 
to  the  pattern  data  in  0.1  increments,  the  amplitudes  of  some  of  the  other  coefficients 
grow  indicating  different  possible  source  locations.  With  a shift  of  0.5  meters  (second  set 
of  coefficients  in  table  2),  three  possible  source  locations  are  indicated  at  -0.5,  + 0.5  and 
+ 1 .5  meters  along  the  X axis.  As  the  shift  is  applied  in  small  steps  (0.  IX  or  less),  the 
possible  source  locations  appear  to  move  around.  When  the  shift  is  sufficient  to  encompass 
the  real  source  locations,  the  amplitudes  predicted  from  the  Fourier  coefficients  reach  a 
maximum.  In  this  case,  an  applied  shift  of  1 .5  meters  (third  set  of  coefficients)  correctly 
predicts  the  source  locations  at  ± 1.5  meters.  If  the  incremental  shifting  is  continued,  the 
amplitudes  will  peak  every  time  they  coincide  with  the  true  source  locations.  For  example, 
when  the  shift  is  2 meters  (last  set  of  coefficients),  the  true  source  locations  are  correctly 
given  by  the  peak  amplitudes. 

The  Fourier  series  expansion  assumes  an  array  of  sources  equally  spaced  at  a half- 
wavelength. It  is  therefore  possible  that  the  array  sources  may  not  exactly  line  up  with 
the  true  source  locations  at  each  increment  of  the  applied  shift,  even  though  the  shift  may 
encompass  the  entire  source  distribution.  This  condition  is  illustrated  in  table  2 for  a shift 
of  1.6  meters  (fourth  set  of  coefficients).  A peak  in  amplitudes  still  occurs  near  the  actual 
source  locations  (±  1.5  meters),  but  other  coefficients  have  significant  amplitudes  as  well. 
Notwithstanding,  the  maximum  amplitudes  occur  when  the  coefficients  align  with  true 
source  locations  and,  again,  they  are  clearly  defined  by  a significant  difference  in  amplitude 
compared  to  the  other  members  of  the  set. 

The  next  example  to  consider  is  that  of  the  array  shown  in  figure  2.  In  this  case, 
the  sources  are  all  located  along  the  positive  X axis.  Figure  2 shows  three  half-wave  dipoles 
located  at  X = 0.5X,  X,  and  4X  in  free  space.  All  three  are  center-fed  with  the  dipole  at 
X = X driven  180°  out  of  phase.  The  pattern  in  the  X-Y  plane  was  calculated  using  NEC 
and  processed  using  the  Prony  and  Fourier  options  of  ADPAP.  These  results  are  given  in 
tables  3 and  4. 

Table  3 gives  the  Prony  residues  and  poles  for  the  cases  of  20.  30,  and  40  poles. 

The  number  of  samples  in  pattern  data  are  consequently  40.  60,  and  80  respectively.  For 
20  poles  (the  first  set  of  poles  in  table  3).  four  possible  source  locations  (pole  numbers  1, 

3,  6,  and  1 1 ) are  identified  by  their  relatively  large  residues.  These  poles  are  located  on 
the  imaginary  axis  (small  real  parts)  of  the  complex  spatial  plane  corresponding  to  the 
correct  source  locations  along  the  array  axis  (see  the  imaginary  part  of  the  poles).  A 
double  pole  is  located  at  3.9816  meters  (poles  6 and  11).  By  increasing  the  number  of 
poles  to  30  (see  second  set  in  table  3),  the  double-pole  ambiguity  is  resolved  and  the 
correct  source  locations  indicated  by  poles  1 . 2.  and  1 7.  Increasing  the  number  of  poles 
to  40  (third  set  in  table  3)  reintroduces  the  double  poles  (pole  numbers  7 and  22).  The 
double  poles  result  from  curve  fitting. 


Table  3.  Prony  residues  and  poles  for  the  pattern  of  the  dipoles  in  figure  2. 


COMPLEX  PRONY 


f REOUENCY 

r 

.2998*009  HZ 

WAVELENGTH 

= 

1.0000  METERS 

NO.  OF 

SAMPLES 

= 

40 

NO.  OF 

POLES  = 

20 

THE  PHASE 

CENTER  IS  LOCATED 

AT  ( 

.0000,  .0000) 

METERS 

RE  SIDUES 

POLES 

MAG 

PHASE 

real 

IMAG 

1 

. 61641*000 

90 

.543 

.1  6647-004 
.16835*000 

.49738*000 

2 

.1 3337-002 

89 

.651 

-.89905*000 

3 

. 84378*000 

144 

.120 

-.22716-002 

.99842*000 

4 

.57368-002 

106 

.658 

.13713*000 

-.29011*001 

5 

.29082-001 

140 

.166 

-.22141-001 

-.46901*001 

* 

-12 

.105 

.507 

.74874-001 

.18799-001 

.39816*001 

-.74497*001 

8 

.40087-002 

-173 

.721 

-.18044*000 

-.16822*001 

9 

.5  1792-002 

159 

.586 

-.45964-001 

.25687*001 

10 

.22496-007 

-167 

.040 

.24731*001 

-.61194*001 

1 1 

.34459*000 

44 

.973 

-.69454-001 

•2^8 15*001 
-.37497*001 

12 

.26857-002 

.93247-003 

-147 

.478 

-.97249-001 

13 

126 

.637 

.57389*000 

.84014*001 

14 

.21924-004 

-130 

.260 

.11939*001 

.62213*001 

15 

.58713-003 

19 

.881 

-.55092*000 

.83520*001 

16 

.37379-001 

-51 

.529 

.52114-002 

-.58396*001 

17 

.51041-004 

-126 

.352 

-.98006*000 

.62740*001 

18 

.4  51  16-002 

5 

.310 

.23334*000 

-.933194001 

19 

.28459-002 

.23647-015 

"if1 

:??  1 

-.17880*000 

-.89795*001 

20 

124 

-.51741*001 

-.47773*001 

r~  ■ - 


Table  3.  Continued. 


NO.  OF  SAMPLES  = 
NO.  OF  POLES  = 


PE  SIDUES 
MAG 


•63641*000 
.841 83*000 
.1 75 43-005 
.1  1646-004 
.85939-003 
.39956-001 
.59911-003 
• 5 8664-004 
.76701-006 
.33275-004 
.1 9747-003 
.92178-002 
.1 5404-006 
.26894-016 
.40842-005 
.1 8212*006 
.664134000 
.4  64  1 9-»002 
.3  28215^004 
.1 2800-005 
.57901-003 
.28726-005 
.97818-003 
.39791-005 

:1K?S:88? 

.4  6884-004 
.49104-005 

:ii*!i:88i 


PHASE 

87.565 

144.866 

84.814 


143.269 
124.645 
164.578 
82.760 
54.174 
140.899 
102.436 
76.108 
166.614 
-86.098 
58.044 
137.201 
28.971 
.v  5 3.44  8 
91  .927 
-12.257 
1Q2.545 
•45.008 
101  .51 3 
58.299 
149.681 
-66.381 
■123.594 
■128.961 


POLES 

PEAL 

.35830-002 

.41111-002 

.20555*001 


004  17.491  .16739*001 
003  -99.899  .85461*000 
001  -146.791  -.11567-001 


.91706*000 
.13957*001 
.20651*001 
.14895*001 
.1 1161*001 
.22092*000 
.20092*001 
.62029*001 
.17667*001 
.18926*001 
-.20902-002 
.39599*000 
.1 3939*001 
.20472*001 
.79053*000 
.17059*001 
.40408*000 
.18667*001 
.94000*000 
.18594*001 
.11360*001 
.16237*001 
.10594*001 
.10580*001 


*..*♦»♦»*»***♦ 


.50659*nnQ 

.99975*000 

-.34274*001 

-.13498*001 

.86623*001 
.10589*002 
-.55086*001 
.12552*002 
.14562*002 
.67071*001 
-.99438*001 
-.18997*001 
-.11880  <002 
-.78345*001 
.40127*001 
-.13776*002 
-.15850*001 
-.50045*001 
.26103*001 
-.90084*001 
.75901*001 
-.32972*001 
.96369*001 
-.67752*001 
.11591*002 
-.10944*002 
-.12866*002 
.13677*002 


Table  3.  Continued. 


NO , or  CAvFLtS  - 
NO.  OF  c 0 L F S - 60 


F £ 5 1 D OF  5 

POLES 

VA0 

PM  A c £ 

PEAL 

IMAG 

1 

.5  5 9 1 7 ♦ 0 0 Q 

*3.165 

_ 

.94046-003 

.47421+000 

«* 

£ 

.46 751-00? 

in v. 704 

.5 1 609+000 

.45050+000 

3 

.*  7647+nOQ 

140.456 

- 

.11*39-001 

.99447+000 

4 

.1 675i-C0? 

-11 6.V5  7 

- 

.434 22-002 

-. 1 ns  1 9 + 001 

5 

.421 55-n0? 

5 7.71V 

.19457-001 

-.<♦4592  + 001 

ft 

.1  8*2*-002 

-1 35. j53 

.65352-001 

-.74896+001 

7 

.3  2231  ♦nrjn 

74.764 

.61469-001 

.39837+001 

a 

.44M3-00? 

-11 3.053 

.4006*-001 

-.67199+001 

9 

,o?1 ^3-004 

-*7.995 

. 75367+000 

-.9*276+001 

10 

.27721-002 

-122.170 

- 

.57737-001 

-.14606+002 

1 1 

.2  g23O-007 

54.97? 

- 

.39743-001 

-.17033+002 

1c 

.2 f 7 V1-OQ3 

-123.436 

.45320+000 

.60126+001 

1 3 

.1  f 0 55-0  32 

60.931 

.74*4  3—001 

-.1*077+002 

14 

. 5 27  90- no  3 

144.444 

.54654+000 

.8031  * +001 

1 5 

.1 C264-O0? 

-123.730 

- 

.1*790-001 

-.331*2+001 

16 

.1 9°  24-00  2 

-127.441 

- 

.29642-001 

.93306+001 

1 7 

.26068-00? 

'6.871 

.6  7 ' 4 5 -00 1 

.10161+002 

1 * 

.ft  (.731-002 

53.562 

- 

.157-4-001 

.124*0+002 

* V 

.a  ?(";0-n32 

170.132 

- 

,5024*  + 000 

. 7*226+000 

JU 

.1  (.n  *«-0C2 

5 7.410 

- 

.26144-001 

-.22013+001 

21 

,oca54-00? 

-1  2 3 . 4 0 1 

.17746-001 

.13943+002 

p 3 

.3>77n»000 

44.394 

- 

.614-1 -001 

.39852+001 

23 

.8^755-003 

-4  3.07  6 

- 

.24952-001 

-.91836+001 

2 4 

,26=41-no? 

- 1 n4 . 2 2 5 

.43771 -002 

.11 326+002 

25 

.1  1214-no? 

40.290 

.73469-001 

-.10859+002 

2 6 

,7r«;46-0C7 

-67.525 

- 

.49794.0OO 

,7*97G+001 

27 

.6  2460-002 

-11  a. 050 

.26747-001 

-.59463+001 

?s 

.36150-002 

-170.731 

.32*71-001 

.14765+002 

? 9 

.31071-002 

-119.665 

,2«9?4-001 

.15645+002 

7 J 

.ID7  50-002 

-122.040 

- 

.35737-002 

-.12093+002 

71 

.2  4Oi,2-00! 

- 1 ■ 6 . 9 9 7 

.1  767  1 +C00 

.1 7520+0n2 

7 2 

.54'<J7-007 

-132.630 

.2O?4'+O01 

. 1 5493^002 

T ■» 

J 

.1  *427-033 

-104.607 

- 

.7  2*  3 0 + 000 

.6n203+001 

u 

.1  5VIJ9-032 

-127.215 

- 

.1  5 *46  + 000 

-.1 '51 2+002 

35 

.1  27<.7-00? 

5 2.65  0 

.10211 -001 

.1*7*4+002 

7 6 

.1  C365-O0v 

7.45  3 

- 

.30*27+001 

-.97130+001 

7 7 

.1  1*63-0 J 2 

-176.373 

- 

.24977-001 

-.13969+302 

7a 

.1  4741-O07 

- 1 7 a . 3 0 0 

• 

.5*786+000 

-.1*566+002 

7 9 

.1041 7-002 

*5.033 

.13504+000 

-.19393+002 

60 

.1  *1  70-00* 

105.363 

“ 

.27291+001 

.17265+002 

Table  4.  Fourier  coefficients  for  the  pattern  of  the  dipoles  in  figure  2. 


FOURIER  COEFFICIENTS  FOR  R E AL < E > 


FREQUENCY  = 

2998*009  HZ 

WAVELENGTH  * 

1.0000  METERS 

NO  . 

OF  SAMPLES  = 

1UU 

NO  . 

OF  TERMS  = 

10 

THE 

PHASE  CFNTER 

IS  LOCATED  AT 

( .0000, 

.0000) 

meters 

N 

scheters) 

I <N> 

J ( N) 

AMP(N  ) 

PHASE 

0 

.0000 

-.32179-002 

-.61  738*000 

.5  71  29-002 

124.282 

1 

.5000 

.22666-001 

.61780*000 

87.897 

2 

3 

1 .0000 

1 .5000 

-. 68859*000 
.65925-002 

-.49685*000 
.1801  2-002 

.84912*000 

.68341-002 

144.188 

-15.281 

4 

2.0000 

-.6  2107-002 

-.14291-002 

.63730-002 

167.041 

5 

2.5000 

.55753-002 

.65068-003 

.561  31-002 

-6.657 

? 

3:9888 

-:15f2i:88f 

-:2??!T-88i 

!63012-OoI 

8 

4.0000 

.35906*000 

-.59092*000 

.69145*000 

58.716 

9 

4.5000 

.89243-003 

.74859-002 

.75390-002 

-83.202 

10 

5.0000 

-.55128-002 

-.36989-002 

.66387-002 

146.139 

! 


Table  4 lists  the  Fourier  coefficients  for  the  pattern  of  the  array  in  figure  2.  The 
dipoles  are  all  nicely  arranged  on  one  side  of  the  phase  center  at  integer  multiples  of  A/2. 

Consequently,  only  one  iteration  of  the  Fourier  series  is  really  necessary  to  identify  the 
sources.  Coefficients  1,  2,  and  8 in  table  4 have  the  largest  magnitudes  and  are  sharply 
defined.  It  is  interesting  to  note  that  the  amplitude  and  phase  of  each  source  predicted 
by  Fourier  is  essentially  the  same  as  those  predicted  by  Prony  (see  the  second  set  of  poles 
in  table  3). 

The  two  examples  thus  far  have  been  nicely  arranged,  ie,  the  dipoles  have  had 
integer  multiples  of  X/2  spacing.  Figure  3 shows  an  array  in  which  two  of  the  monopoles 
are  OJA  apart.  The  monopoles  are  all  located  on  one  side  of  the  phase  center  as  in  the 
previous  case.  For  this  example,  the  patterns  were  calculated  for  the  case  with  only  one 
monopole  excited:  the  one  at  X = 0.5A.  In  this  case,  the  monopole  at  1 .5X  is  too  weakly 
coupled  to  contribute  much  to  the  pattern.  The  results  of  applying  the  Prony  and  Fourier 
processing  are  given  in  tables  5 and  6. 

Table  5 shows  the  Prony  residues  and  poles  for  the  array  of  monopoles  in  figure 
3.  After  a number  of  iterations  using  a different  number  of  poles  each  time,  it  was  found 
that  about  30  poles  with  60  sample  points  is  satisfactory.  Two  poles  are  distinctly  identified 
by  their  relatively  larger  residues  (pole  numbers  2 and  3).  They  are  located  on  the  imaginary 
axis  of  the  complex  spatial  plane  as  required.  Pole  number  2 has  a greater  residue  than  pole 
number  3,  which  should  be  expected  since  it  corresponds  to  the  driven  source  at  0.5  meters 
on  the  array  axis.  The  Prony  process  has  no  difficulty  in  correctly  resolving  the  source 
locations  even  though  they  are  separated  by  only  0.3A. 

Table  6 shows  four  successive  iterations  of  the  Fourier  process  in  which  the  phase  , I 

center  is  shifted  0. 1 meter  along  the  array  axis  each  time.  The  amplitudes  peak  around  the 
proper  locations  for  two  sources  at  X = 0.5A  and  0.8A.  However,  the  true  locations  are  not 
readily  apparent.  The  Fourier  process  cannot  resolve  sources  that  are  closer  than  0.5A  apart. 

3.2  ANTENNAS  IN  A SIMPLE  ENVIRONMENT 

The  previous  three  examples  were  for  simple  arrays.  The  only  radiating  structures 
were  the  dipoles,  and  they  were  arranged  along  a common  axis. 

Figure  4 shows  an  array  of  35-foot  whips  mounted  on  a box  whose  dimensions  are 
similar  to  the  superstructure  of  a small  ship.  The  35-foot  whips  are  8.69  meters  apart  and 
form  an  array  whose  axis  intercepts  the  X-Z  plane  at  164.75°  (see  figure  5).  At  hf  frequencies 
where  these  antennas  are  employed,  currents,  which  also  radiate  and  influence  the  patterns, 
are  induced  on  the  box  surface.  For  this  illustration,  the  whip  at  the  box  center  is  driven 
at  20  MHz  so  the  whips  are  0.58A  apart.  Both  Prony  and  Fourier  should  be  able  to  resolve 
the  true  sources.  The  whips  and  the  box  were  modeled  using  NEC  and  the  resulting  pattern 
processed  using  ADPAP.  Some  of  these  results  are  listed  in  tables  7 and  8. 

The  patterns  calculated  by  NEC  were  for  the  phase  center  at  the  origin  as  shown  in 
figure  4.  When  Prony  was  applied,  the  resulting  matrices  were  ill-conditioned,  causing  an 
abnormal  abort  of  the  program.  This  problem  was  observed  even  with  the  simple  array 
problems  previously  discussed  when  the  assumed  array  axis  did  not  coincide  with  the  true 
axis  of  the  array.  If  the  assumed  axis  is  off  5 to  6 degrees,  acceptable  results  may  be  obtained. 
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Table  5.  Prony  residues  and  poles  for  the  pattern  of  the  array  in  figure  3. 


COMPLEX  PRONY 


FREQUENCY  = 
WAVELENGTH  = 


.2998*009  HZ 

1.0000  METERS 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 
1 1 
12 

11 

11 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 


NO  . 

OF 

SAMPLES  = 

60 

NO  . 

OF 

POLES  = 

30 

PHASE  CENTER  IS  LOCATED  AT  ( 

.0000.  . 

0000)  METERS 

RESIDUES 

POLES 

MAG 

PHASE 

REAL 

IMA6 

.42136-003 

92.977 

-.35895-001 

-.1  3346*0nn 
.50336*000 

.1 0167*001 

52.877 

-.16557-002 

:12?8U882 

152.947 

.92165-002 

.80206*000 

114.411 

.53449*000 

-.16028*001 

.33125-004 

.10081-005 

-133.740 

-30.790 

.4394  3*nnn 
.1001  2*001 

.35144*001 

.79483*001 

-.43460*001 

.1 6913-004 
.89118-001 

-25.596 

.60108*000 

-122.082 

.19033-002 

.14921*001 

.56485-021 

-90.995 

.69601*001 

-.76495*001 

.1 3287-004 

-36.419 

.58484*000 

.12679*002 

.5  7516-005 

-143.455 

.85356*000 

-.70479*001 

.1 7603-004 

-53.222 

-.76860*000 

-.37654*001 

.1 2042-005 

-5  1.646 

.1  1626*001 

-.91657*001 

.40049-004 

14.728 

.26872*000 

.59783*001 

.48042-004 

88.848 

-.32948*000 

-.11095*002 

.44766-005 

-121 .928 

.75230*000 

-.12900*002 

.57594-004 

-100.471 

-.49487*000 

-.20345*001 

.50265-004 
.1 4092-004 
.79368-005 

29.142 

-.50913*000 

.39355*001 

33.361 

-.82421*000 

.77080*001 

-3  0.99  9 

-.86774*000 

-.72533*001 

.64959-003 

-32.634 

-.12842*000 

.97510*001 

.1 9237-003 

-158.660 

-.45028*000 

.99775*001 

.1 3309-004 

-42.618 

-.85969*000 

-.55038*001 

.38634-004 

.65818-005 

16.165 

48.248 

-.58138*000 

-.77152*000 

.58294*001 

.11343*002 

.14786-004 

-72.243 

-.678 32*000 

.14463*002 

.1 5759-003 

-35.370 

.91925-001 

.14376*002 

2.161 

90.628 

-.11480*001 

-.75441*000 

-.91085*001 

-.13308*002 

.1 8148-01 1 

1 1.255 

-.34657*001 

-.11  840*002 
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Table  6.  Fourier  coefficients  for  the  pattern  of  the  array  in  figure  3. 


FOURIER  COEFFICIENTS  FOR  REAL(E) 


FREOUENCY  = 

. 2998*009  H2 

WAVELENGTH  = 

1.0000  METERS 

NO.  OF  SAMPLES 

= 100 

NO.  OF  TERMS  = 

10 

THE 

PHASE  CENTER 

HAS  BEEN  SHIFTED 

BY  .0000 

METERS  along 

ARRAY  AXIS 

N 

scmeters) 

I (N) 

J (N) 

AM P ( N ) 

PH  ASE 

0 

.0Q00 

.76633-001 

.93782-001 

-35  .200 

l 

i :5888 

::i?5$$:888 

:!5JIS:888 

3 

1 .5000 

.34032-001 

.1441 6*000 

.14812*000 

-76.717 

4 

2.0000 

-.46162-001 

-.45243-001 

.64636-001 

135.576 

5 

2.5000 

.33406-001 

.33733-001 

.47475-001 

-45  .279 

6 

3.0000 

-.27106-001 

- .27066-001 

.38306-001 

1 35  .043 

7 

3.5000 

.23517-001 

.22632-001 

.32638-001 

-43  .902 

8 

4.0000 

-.21260-001 

-.19483-001 

.28838-001 

137.497 

9 

4.5000 

.19725-001 

.1  7085-001 

.26095-001 

.24005-001 

-40.898 

10 

5.0000 

-.18637-001 

-.15131-001 

140.928 

*********** 

* * * * 

THE 

PHASE  CENTER 

HAS  BEEN  SHIFTED 

BY  .1000 

METERS  ALONG 

array  AXIS 

N 

SCMETERS) 

I (N) 

J ( N) 

AMP  ( N ) 

phase 

0 

- . 1 Q00 

-.54365-001 

.1 5809*000 

-110.113 

1 

.4000 

.49828*000 

-.70895*000 

.86654*000 

54  .899 

2 

.9000 

-.29279*000 

- .4861  3*000 

.56749*000 

121 .060 

3 

1 .4000 

-.13045-001 

.2321 5*000 

.23252*000 

.96154-001 

-93  .216 

4 

1 .9000 

-.30836-001 

-.91  075-001 

108  .705 

5 

6 

1:1888 

-:HM:881 

.7431 4-001 
-.61  337-001 

-73.197 

107.958 

7 

1:1888 

.18710-001 

-.18062-001 

.5201 8-001 

.55280-001 

-70.217 

8 

-.45  100-001 

.48582-001 

111  .825 

9 

4.4000 

.17641-001 

.39726-001 

.4  3467-001 

-66.055 

10 

4.9000 

-.17357-001 

-.35368-001 

.39397-001 

116.139 

Table  6.  Continued. 


THE  PHASE  CENTER  HAS  BEEN  SHITTED  BY  .2000  METERS  AL0N6  ARRAY  AXIS 


N 

S (METE  RS  ) 

I (N) 

J ( N ) 

AMP  ( N ) 

PHASE 

0 

.8000 

-.13361*000 

.21112*000 

-129.261 

1 

2 

,5405l*nnn 
-.18921 *000 

- .49263*000 

- .7302  3 + 0nn 
.25506*000 

.73132*000 

.75435*000 

42.346 

104.526 

3 

1 .3000 

-.98411-001 

.27340*000 

-111 .097 

4 

1.8000 

2.3000 

.26219-00? 

-.95415-001 

.95451-001 

88.426 

5 

-.27821-002 

.86249-001 

.86294-001 

-91  .848 

6 

7 

2.8000 

3.3000 

-.11809-002 

.39743-002 

-.7231  9-001 
.61686-001 

.72329-001 

.61813-001 

90.936 

-86.314 

8 

3.8000 

-.58690-002 

-.5  362  1-001 

.53941-00 1 

96.246 

9 

4.3000 

.71759-002 

.47300-001 

.4  7841-001 

-81 .373 

10 

4.8000 

-.81162-002 

- .42167-001 

.42941-001 

100.895 

* 


THE 

PHASE  CENTER 

HAS  BEEN  SHIFTED 

BY  .3000 

METERS  AL0N6 

ARRAY  AXIS 

N 

S(METERS) 

I (N) 

J ( N) 

A"P( N) 

PHASE 

0 

-.3000 

-.15062*000 

.1  9083  + 000 

-142.120 

1 

.2000 

.46774*000 

-.27750*000 

.54387+000 

30.680 

2 

.7000 

-.16414-001 

-.89593*000 

.89608*000 

91  .050 

3 

1 .2000 
1.7000 

-.20120*000 

.18792*000 

.27531*000 

-136.954 

4 

.39401-001 

- .49406-001 

.63194-001 

51  .428 

5 

6 

i:?888 

-.3  23Q7-QQ1 
.21383-001 

-iHSJUBBI 

.71563-001 

.59200-001 

-116.837 
68  .826 

7 

3.2000 

-.14630-001 

.47559-001 

.4  9759-001 

-107.099 

8 

3.7000 

.10290-001 

-.41  534-001 

.42789-001 

76.085 

9 

4.2000 

-.73599-002 

.36733-001 

.37463-001 

-101 .330 

10 

4.7000 

.52885-002 

-.32798-001 

.33222-001 

80  .840 

I 
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Table  7.  Prony  residues  and  poles  for  the  pattern  produced  by  the  configuration  in  figure  4. 

COMPLEX  P 3 ON  Y 


FREQUENCY  = 

. 200?+00p  h 2 

WAVELENGTH  = 

1 A .9900  VE  T F»  S 

START  = 165.000 

DEGREES 

NO.  OF  SAMPLES  = 

5 0 

NO.  OF  POLES  = 

2 5 

THE 

PHASE  CENTER  IS  LOCATED  AT  ( 

9.1440,  ’. 

0430)  METERS 

RESIDUES 

w A 0 

PHASE 

POLES 

real 

IMAG 

1 

.2  5221  + 002 

-1 6.960 

-.53407+000 

.1  * 749-nn  1 
-.64  258  + 0(31 

2 

•2*126+001 

-2  9.669 

-.10920+001 

3 

.«  ?63O»001 

30.663 

-.29286+000 

.54185+001 

A 

.5  5 ’ 81 -00  3 

13.573 

.67353+001 

-.26556+002 

5 

.5  5n  54-00  2 

-106.777 

.4371 0+001 
.83893+001 

.67715+002 

6 

.3 5237-002 

6 0. 9 « 4 

-.11147+003 

-.14726+003 

7 

.3t5:s5-003 

-p  9 . 64  2 

.99448+001 

6 

.875b2-03T 

-77  .<.46 

.1  7931  +0C3 

.20237+002 

9 

.1  ’759-noi 

47.302 

-.26526+001 

.1*104+003 

10 

.4  tl  1 2-002 

-4  7.599 

.46636+001 

.13505+003 

1 1 

.2 1204-001 

-R5.734 

-.16040+001 

-.79451+002 

1 2 

.*7169-006 

113.407 

.34350+002 

-.9*573+001 

1 7 

.7  6*27-003 

-173.684 

.1 1648*001 

.37464+002 

1 4 

. 2 5227-001 

-’4.403 

ziiiuizm 

.77511*00 2 

.10205+003 

1 5 

.4 c7 11-001 

’4.755 

1 6 

.1 C4i^_poi 

-1*1.216 

. 50439*000 

.11738+003 

1 7 

.71247-002 

-55.854 

-.16*85+001 

.15175+003 

1 6 

.1  <520-001 

29.833 

-.1 0278+C01 

-.16403+003 

1 9 

.4  ’051-003 

125.540 

-.1 2*33+002 

-.14635+003 

*'  r\ 

L.  Ns 

.33737-001 

33.516 

-.70491+000 

-.70479+002 

21 

.72450-001 

’3.648 

-.84397*000 

-.81545+002 

22 

.3  43  35-001 

’3.810 

-.24770+001 

.91 324+002 

2 3 

.1 7263-005 

-60.452 

-.’3251+002 

-.64318+001 

24 

.27559-003 

-*7.137 

-.16233+002 

-.11272+003 

25 

.2(471-002 

-114.743 

-.90460+001 

.16087+003 

1 


[q 


Table  8.  Fourier  coefficients  for  the  pattern  of  the  configuration  in  rigure  4. 


FOURIER  COEFFICIENTS  F OR  REAL<E) 


FREQUENCY  = .2000*008  HZ 

WAVELENGTH  = 14, 

.9900  *ETERS 

START 

= .DUO  DEGREES 

NO  . OF 
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Larger  angles  often  cause  abnormal  program  termination.  Another  complication  is  that  for 
the  problem  of  figure  4,  the  sources  are  not  all  on  a single  axis  since  the  box  itself  also 
radiates.  Misalignment  of  multiple  sources  can  be  thought  of  as  noise  in  the  pattern  data 
and  the  Prony  process  is  known  to  be  very  susceptible  to  noise. 

The  assumed  array  axis  is  determined  by  the  initial  angle  in  the  pattern  data.  If  the 
pattern  data  is  in  the  normal  order,  the  initial  angle  is  zero  degrees  and  the  assumed  array 
axis  lies  in  this  direction,  ie,  along  the  X axis.  If  the  pattern  data  were  re-ordered  so  that 
the  initial  angle  is  90  degrees,  for  example,  the  assumed  array  axis  would  coincide  with  that 
angle  - the  Y axis.  Thus  in  ADPAP,  the  assumed  array  axis  may  be  changed  by  selecting 
the  starting  angle  of  the  pattern  to  be  processed.  The  angle  specified  must  coincide  with 
an  angle  in  the  pattern  data;  eg,  if  the  pattern  is  specified  every  two  degrees  starting  at 
zero,  only  even  angles  may  be  used  to  orient  the  assumed  array  axis. 

Making  use  of  the  axis  orientation  option  as  well  as  the  shifting  theorem  option  ot 
ADPAP,  the  Prony  residues  and  poles  listed  in  table  7 may  be  obtained.  The  array  axis  was 
rotated  165  degrees  to  align  the  assumed  array  axis  to  the  nearest  whole  angle  of  the  true 
array  axis  shown  in  figure  5 (the  pattern  data  was  calculated  by  NEC  in  5 increments). 

The  phase  center  has  also  been  shifted  to  coincide  with  the  location  of  the  35-foot  whip 
at  the  box  center. 

Only  the  first  three  poles  in  table  7 have  significant  residues  with  the  number  1 and 
3 poles  close  to  the  imaginary  axis.  Thus  two  sources  are  identified:  one  at  the  center  of 
the  assumed  anay  axis  and  the  other  at  S = +8.4 1 85  meters.  Further  iterations  of  Prony 
with  more  poles  lead  to  the  same  results.  Noise  introduced  by  misalignment  of  the  array 
axis  and  radiation  from  parts  of  the  box  are  responsible  for  the  large  real  parts  of  the  poles. 
Most  of  the  errors  in  pole  location  and  source  strength  are  attributed  to  the  same  cause.  When 
the  problem  is  over  specified  by  requiring  too  many  poles,  additional  curve  fitting  sources 
are  introduced  which  are  not  clearly  different  from  the  true  source.  Thus  a possible  source 
is  indicated  by  pole  number  2 at  S = -8.4258  meters. 

Table  8 lists  Fourier  coefficients  for  the  pattern  of  the  configuration  in  figure  4. 

The  first  set  of  coefficients  is  for  direct  application  of  the  Fourier  process  to  the  pattern 
data  as  it  was  originally  generated  (the  phase  center  is  at  the  origin  in  figure  4).  Multiple 
sources  are  indicated  between  the  origin  and  29.98  meters  along  the  X axis.  Iteration  of 
the  shifting  theorem  does  not  clearly  identify  the  sources  when  the  assumed  array  axis  is 
badly  misaligned.  The  second  set  of  coefficients  in  table  8 results  from  reorienting  the 
assumed  array  axis  to  165°  and  shifting  the  phase  center  to  the  center  of  the  box,  as  was 
the  case  for  Prony.  After  incrementally  shifting  along  the  array  axis,  the  maximum  ampli- 
tudes were  found  with  a shift  of  8.69  meters.  The  two  sources  are  indicated  by  coefficients 
1 and  2 at  a halt-wavelength  separation.  The  strongest  source  is  at  the  array  center,  which 
is  the  location  of  the  driven  antenna.  Again  there  is  good  correlation  between  Fourier  and 
Prony  for  the  source  amplitude  and  phase. 

3.3  ANTENNAS  ON  A SHIP 

Figure  6 shows  a sketch  of  a wire  grid  model  of  the  PGG,  a small  ship  about  58 
meters  long.  The  two  35-foot  whips  are  about  23.24  meters  apart.  Figure  7 shows  a top 
view  of  the  PGG  giving  exact  antenna  locations. 


NEC  was  used  to  model  the  35-foot  whips  and  the  ship,  using  a wire  grid  scheme 
for  conducting  surfaces  as  shown  in  figure  6.  Patterns  were  calculated  using  the  origin 
shown  as  the  phase  center.  Azimuth  patterns  were  determined  in  5-degree  increments 
for  a 5-degree  elevation  angle.  The  pattern  of  the  aft-most  whip  (antenna  2-1 ) at  10  MHz 
was  processed  using  ADPAP.  At  10  MHz,  the  antennas  are  sufficiently  separated  for 
resolution  by  both  Prony  and  Fourier  processing.  Some  of  these  results  are  shown  in 
tables  9 and  10. 

Table  9 lists  the  Prony  residues  and  poles  for  the  10-MHz  pattern  of  antenna  2-1. 
The  assumed  array  axis  is  set  at  5 degrees  to  coincide,  within  the  nearest  whole  degree  of 
pattern  data,  to  the  actual  alignment  of  the  whips  (see  figure  7).  The  phase  center  was 
not  moved  from  the  origin  defined  in  figures  6 and  7.  A number  of  iterations  of  the 
Prony  process  for  10  through  50  poles  were  performed.  The  results  of  table  9 are  typical. 
Two  poles  (pole  numbers  3 and  4)  have  relatively  large  residues  and  are  located  at  13.3 
and  36.5  meters  along  the  assumed  array  axis.  The  whip  separation  is  correctly  predicted. 
A small  error  in  specific  location  is  seen.  Also,  the  poles  are  not  as  close  to  the  imaginary 
axis  in  the  complex  spatial  domain  as  they  were  for  the  simpler  arrays.  Both  discrepancies 
arc  attributable  to  noise  introduced  by  (1)  misalignment  of  the  assumed  array  axis  with 
the  true  sources  and  (2)  radiation  from  structural  members  of  the  ship. 

Table  10  lists  the  Fourier  coefficients  derived  from  the  pattern  of  antenna  2-1. 

The  assumed  array  axis  is  set  at  5 degrees,  as  in  the  case  of  Prony  processing.  The  shifting 
option  is  applied  in  small  increments  to  shift  the  phase  center  along  the  array  axis  in  search 
of  peaks  in  the  amplitudes.  Three  of  the  iterations  are  shown  to  illustrate  this  approach. 
One  peak  occurs  at  1 1.99  meters  (coefficient  number  1 in  set  2)  and  another  at  23.994 
meters  (coefficient  number  2 in  set  3).  The  locations  (S)  are  with  respect  to  the  original 
phase  center.  If  smaller  increments  in  the  shifting  theorem  are  used,  a more  exact  deter- 
mination of  the  source  locations  is  possible.  This  process  is  tedious  and  may  require  an 
excessive  number  of  iterations  to  equal  the  accuracy  of  Prony  processing. 

3.4  SUMMARY 

The  preceding  sections  have  described  how  Prony  processing  may  be  applied  to 
antenna  patterns  to  determine  source  locations.  Section  3.3  demonstrated  that  Prony 
processing  can  accurately  identify  the  principal  radiators  even  for  the  cluttered  shipboard 
environment.  This  technique  could  be  a very  useful  tool  during  the  shipboard  antenna 
arrangement  design  in  cases  where  one  or  more  parasitic  structures  are  present.  Not  only 
can  these  structures  be  located  by  the  Prony  poles,  but  their  relative  importance  can  be 
deduced  from  the  size  of  the  residues. 

Prony  processing  is  not  the  only  possible  approach,  as  demonstrated  via  Fourier 
series  processing.  But  the  chief  advantages  of  Prony  are  that  ( 1 ) unlike  the  Fourier  series 
approach,  Prony  processing  can  resolve  sources  separated  by  less  than  one-half  wavelength, 
and  (2)  the  solution  is  found  with  fewer  iterations.  Perhaps  with  experience  or  some 
suitable  eigenvalue  technique  (suggested  by  Lawrence  Livermore  Laboratory)^,  only  one 
or  two  iterations  may  be  required.  The  only  disadvantage  of  Prony  is  its  susceptibility 
to  noise  in  the  pattern  data.  This  problem  is  now  receiving  considerable  attention  from 
the  technical  community. 
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Tabic  9.  Prony  residues  and  poles  for  the  pattern  of  the  configuration  in  figure  6. 
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Table  10.  Fourier  coefficients  for  the  pattern  of  the  configuration  of  figure  6. 


FOURIER  COEFFICIENTS  FOR  REAL(E) 

FREQUENCY  = .1COO*0OP  HZ 

•AVFLENGTH  = 20.0800  ME  T E R S 

STAPT  = 5.000  DE6RFES 

NO.  OF  SAMPLES  = 100 

NO.  OF  TERMS  = 10 


THF  PHASE  CENTER  IS  LOCATED  AT  ( .0000,  .0000)  METERS 


N 

SC  METERS) 

I (N) 

J ( N) 

AMp (N  ) 

PH  ASE 

0 

.0090 

-.30765-001 

.46274+000 

-93.812 

1 

1 4 . 9900 

.10464+001 

- .31 31 5+000 

.1  0922  + 001 

16.661 

2 

2 0. op  no 

-.74211+009 

- .1 735 f ♦000 

.38363+000 

153.09? 

3 

4 4 . 97  nn 

.42590-001 

.2  763  2 + 000 

.24062*000 

-79.805 

4 

59.9600 

-.56157-301 

- .1  5 45  5 + 000 

.1 6444+000 

109.969 

5 

74.9500 

.49434-001 

.1  293  2 + 000 

.1 3845+000 

-69.081 

0 

: 9 . 94  00 

-.45752-001 

- .1  036  5 + 000 

.11330+000 

113.818 

7 

134.9300 

.42355-001 

.3  6 33  4-001 

.96164-001 

-63.867 

St 

119.9200 

-.79845-001 

- .73972-001 

.83986-001 

118 .322 

9 

1 3 4 .91  no 

.3*049-031 

.64499-001 

,74s  86-r01 

-59.463 

1 0 

149.9000 

-.76697-301 

- .571  7 6-001 

.6  79  U6- 001 

122.712 

*************** 

the  phase  cfnter  has  pefn  shitted  py  2.09 so 
N c(“ETERS)  I < N ) J ( \ ) 

Q -2. OV80  -.31674+CCO 

1 11.0920  . 1 2864  + 001  -.64409-001 

2 26.0820  -. 24370+nC0  -.42  81  1 + 000 

3 4 1 . 0720  -.  1 9037-001  .3662  3*GC0 

4 56.0620  -.42227-001  -.22204+000 

5 71.0520  .44513-001  .17641+003 

6 :6.O4  20  -.44877-001  - .1  778  2*000 

7 131.9320  .43462-001  .11365+000 

8 11 6 . °2  2C  -.42  261  -001  -. 96656-301 

V 1 3 1 . 91  20  .4  1 76s  -001  .8  4 32  9-001 

10  146.9020  -.40647-031  -.74634-001 


*************** 


THE 

PHASE  CENTER 

HAS  OEEN  SHI ETED 

b Y 5.99  60 

METERS  along 

array  AXIS 

N 

S ( M e T E R S ) 

I ( N ) 

J ( N ) 

AMP(N  ) 

PHASE 

J 

-5.99*0 

-.45619+000 

.52492*000 

-150.351 

1 

w . 9 9 4 0 

. 1 264  2*901 

.1  5430+000 

.1 2736*001 

-6.958 

L, 

27.984 0 

.1 1 777-OOI 

- .6004  1 + 000 

.60052*000 

68.914 

! j 

3 » .9740 

-.14114*000 

.39577+000 

.42018*000 

-109.627 

x 4 

5 1.9640 

.42724-9Q7 

- .2983  0*000 

.20835*000 

68.825 

. 5 

6 = .95  40 

. 1 1987-001 

.1  *00  5*000 

.1 6049+003 

-85.717 

• 6 

'*3,9440 

-.19792-001 

- .12030*000 

.12192+000 

99.343 

7 

v » . 0 i L 0 

.226  57-001 

.9791 *-001 

.1 0050+000 

-76 .971 

*8 

1 1 3 . 92  40 

-.2440“-001 

- .8  295  4-001 

.86470-001 

106.396 

9 

1 : r .91 40 

. 2 5 599-00-1 

.72024-001 

.7&438-C01 

-70.43? 

10 

143.9040 

-.26773-001 

- .6 7 66 t -001 

.68914-001 

112.501 

29 


METERS  ALONG  ARRAY  AXIS 


ampcn  ) 


PHASE 


.5 101»*000 
.1 2*80+001 
.49261+003 
.36677+000 
.22602+000 
.1  0 1 94  + 000 
.1  4494  + 000 
.1 216**000 
.1 0568+000 
.93929-301 
.8493T-001 


•1  28.378 
2 .866 
119.650 
-93.116 
100.768 
-75.53? 
108.037 
-69.073 
113.572 
-63.870 
118.571 


4.0  APPLICATION  TO  PATTERN  DATA  COMPRESSION 

Prony  processing  of  pattern  data  to  obtain  source  location  and  strength  information 
requires  both  magnitude  and  phase  of  the  pattern.  The  pattern  data  used  in  the  examples 
of  section  3.0  consist  of  from  72  to  360  complex  numbers,  depending  on  whether  data  are 
supplied  in  5-  or  1 -degree  increments.  For  acceptable  accuracy,  anywhere  from  20  to  40 
poles  were  determined. 

If,  on  the  other  hand,  the  source  location  and  strength  information  is  not  required, 
but  instead  the  objective  is  to  reproduce  the  original  pattern  data  from  the  Prony  derived 
exponential  series,  fewer  poles  will  be  required.  Either  fewer  poles  need  be  determined  in 
the  first  place  or  some  of  the  poles  of  a set  may  be  truncated  from  the  series.  The  number 
of  poles  required  for  adequate  representation  is  determined  in  part  by  the  application  and 
in  part  by  trial.  In  fact,  the  original  pattern  data  need  not  include  the  phase,  if  this  infor- 
mation is  not  required  by  the  application.  Of  course,  the  Prony  residues  and  poles  would 
lose  their  physical  identities,  but  may  still  be  used  to  retrieve  the  original  pattern  data. 

Large  volumes  of  pattern  data  are  already  measured  on  the  NOSC  pattern  range  as 
a routine  part  of  shipboard  antenna  arrangement  procedures.  (These  data  do  not  as  a rule 
include  phase  information.)  The  pattern  data  are  used  immediately  at  the  time  of  collection 
for  some  statistical  processing,  some  are  plotted  and  all  may  be  stored  tor  future  reference 
on  magnetic  tape.  Representation  of  these  data  as  sets  of  Prony  derived  complex  numbers 
could  significantly  reduce  storage  requirements.  (However,  a Fourier  series  representation 
could  be  equally  as  good  in  this  case.)  Instead  ot  retaining  360  bits  ol  information  per 
pattern,  only  a few  pairs  of  complex  numbers  need  be  stored,  depending  on  the  required 
accuracy. 

Another  application  of  Prony  for  pattern  storage  and  retrieval  is  in  direction  finding 
(DF)  using  adaptive  arrays.  Especially  in  the  case  where  an  adaptive  array  uses  a number  ol 
drastically  dissimilar  antennas,  the  Prony  representation  of  antenna  patterns  could  be  useful. 

An  adaptive  array  automatically  adjusts  the  phase  and  amplitude  of  each  receive 
channel  to  optimize  the  composite  signal  to  some  predetermined  criterion,  usually  signal 
to-noise  ratio.  The  array  processor  makes  no  use  of  antenna  locations,  antenna  patterns, 
and  phase  delays  in  coaxial  lines  and  multicoupler  filters.  Hence,  the  processor  generally 
does  not  have  enough  information  to  determine  the  array  pattern.  For  purposes  ol  adaptive 
nulling  out  of  interfering  signals,  the  exact  antenna  pattern  shape  is  unnecessary  information, 
but  it  is  essential  for  a DF  capability. 

The  information  available  from  the  array  processor  after  it  has  adapted  is  a set  ot 
complex  weights  (amplitude  and  phase)  for  each  signal  channel.  This  information  could 
be  combined  with  the  phase  delay  in  each  channel  path  plus  accurate  antenna  patterns,  to 
be  derived  from  Prony  poles  and  residues,  and  interpreted  for  DF  information.  Based  on 
a recent  paper  by  Dr.  Kazuki  Takao  and  Dr.  S.  Keiichi.*0  it  may  also  be  possible  to 
distinguish  between  ground  and  skywave  signals  or  other  multipath  signals  using  Prony 
processing. 

10.  Faculty  of  Engineering.  Kyoto  University,  Kyoto.  Japan.  EMCJ  77-34,  Measurement  ot  TV  Ghost 
Waves  by  a Small  Antenna,  by  K.  Takao  and  S.  Keiichi,  28  October  1477 


Another  useful  application  of  Prony  processing  is  to  couple  two  numerical  antenna 
modeling  techniques,  method  of  moments  and  the  geometric  theory  of  diffraction  (GTD). 

A general  purpose  electromagnetics  code  for  communications  antennas  which  employs  GTD 
is  under  development  at  Ohio  State  University  (OSU).1 1 This  code  computes  the  fields  of 
one  or  more  dipoles  including  the  effects  of  nearby  scattering  objects  composed  of  flat 
plates  and  cylinders.  The  code  is  valid  at  frequencies  where  asymptotic  techniques  are 
appropriate.  NEC,  which  has  previously  been  described,  is  a general  purpose  antenna  code 
which  solves  the  boundary  value  problem  for  the  antenna  currents  using  the  method  of 
moments.  The  cost  of  solving  for  the  currents  on  an  antenna  and  all  appropriate  structures 
using  NEC  soon  becomes  unreasonable  as  the  frequency  is  increased.  The  solution  is  to 
switch  to  an  asymtotic  technique  like  the  GTD  code  of  OSU.  However,  for  some  antennas, 
a simple  dipole  representation  and  an  equivalent  dipole  array  may  not  be  readily  obvious. 
Using  NEC,  the  antenna  can  be  accurately  modeled  by  itself;  the  current  distribution  is 
found  and  the  antenna  pattern  calculated  (in  free  space).  Prony  processing  of  the  pattern 
will  then  replace  the  antenna  currents  with  an  array  of  dipoles  whose  strengths  (both  mag- 
nitude and  phase)  and  locations  are  specified  by  the  residues  and  poles.  These  sources  are 
then  used  in  the  GTD  code  to  calculate  the  patterns,  including  the  scattering  from  nearby 
structures.  This  technique  has  shown  a considerable  savings  in  computer  costs  for  two  test 
cases  (not  shown  here)  and  warrants  further  study. 

5.0  CONCLUSIONS 

The  technical  objective  of  this  study  is  to  assess  the  feasibility  and  the  utility  of 
inverting  the  far  field  patterns  of  shipboard  antennas.  The  feasibility  of  using  Prony’s 
method  is  illustrated  in  sections  2.0  and  3.0.  Section  2.0  describes  the  Prony  method 
for  deriving  the  coefficients  and  arguments  of  an  exponential  series  representation  of 
far  field  antenna  patterns.  These  coefficients  and  arguments  can  be  interpreted  to  reveal 
source  strength  and  location  as  demonstrated  in  section  3.0.  The  question  of  utility  of 
the  Prony  transformation  or  inversion  is  not  as  easily  answered  as  is  the  question  of 
feasibility.  Two  applications  were  considered;  data  characterization  and  data  compression. 

An  application  of  Prony’s  method  for  data  characterization  is  the  imaging  of 
radiation  sources.  Section  3.0  demonstrates  that  the  Prony  derived  coefficients  and 
arguments  can  provide  source  strength  and  location  information  for  a given  shipboard 
antenna  pattern.  This  is  possible  because  of  the  exponential  relation  between  an  array 
of  sources  and  their  pattern.  The  usefulness  of  this  approach  for  identifying  shipboard 
sources  from  their  patterns  is  limited  by  two  disadvantages  of  Prony’s  method.  The 
principal  disadvantage  of  Prony’s  method  for  data  characterization  is  its  susceptibility 
to  noise  in  the  pattern  data.  Noise  introduces  errors  in  the  source  representation.  The 
noise  in  the  pattern  data  is  introduced  by  the  misalignment  of  sources  and/or  the  dis- 
tributed nature  of  radiation  from  deliberate  antennas  with  parasites.  The  second  dis- 
advantage is  that  the  interpretation  of  the  Prony  coefficients  and  arguments  is  not  always 
straightforward.  Considerable  experience  with  the  application  of  Prony  processing  and 
some  knowledge  of  the  primary  source  locations  are  prerequisites.  Even  then  the  pro- 
cess of  interpretation  of  the  Prony  data  is  not  always  explicit.  As  a result,  the  cost  of 
the  engineering  time  far  outweighs  the  cost  of  the  computational  time. 


1 1 . Ohio  State  University  Electroseience  Laboratory . TR  4508-8  (784  508),  User's  Manual  for 
Plates  and  Cylinder  Computer  Code,  by  RJ  Marhefka,  March  1978. 


Of  course,  Prony’s  method  is  not  the  sole  technique  for  successfully  extracting 
source  information  from  pattern  data.  One  other  approach,  also  demonstrated  in  section 
3.0,  is  the  use  of  a truncated  Fourier  series.  The  accuracy  of  this  approach  is  not  as 
sensitive  to  noisy  data  and  there  is  no  difference  in  computation  costs.  The  principal 
disadvantage  of  the  Fourier  series  approach  is  that  resolution  of  closely  spaced  sources 
is,  at  best,  one-half  wavelength.  However,  in  one  example  of  section  3.0,  Prony’s  method 
was  able  to  resolve  two  sources  separated  by  only  one-third  wavelength. 

Another  application  of  Prony’s  method  for  data  characterization  is  that  of  develop- 
ing a model  for  a radiating  platform  that,  in  a perturbation  sense,  could  be  used  to  predict 
the  new  behavior  when  the  geometry  is  modified.  For  example,  the  data  from  the  NOSC 
scale  model  range  could  be  used  to  design  and  evaluate  antenna  systems  for  ships,  especially 
when  it  is  not  practical  to  model  all  possible  configurations.  However,  the  noise  problem 
and  interpretation  problem  severely  limit  the  use  of  Prony’s  method  in  shipboard  antenna 
design.  Further  development  will  be  necessary  before  Prony’s  method  becomes  a practical 
tool  for  the  shipboard  antenna  designer.  Since  further  development  of  Prony’s  method  is 
being  pursued  under  other  applications,  it  is  appropriate  to  await  possible  improvements 
from  these  efforts  rather  than  to  continue  in  this  one. 

Prony’s  method  can  also  be  applied  to  data  compression.  When  the  objective  is 
data  storage  and  retrieval  of  large  volumes  of  data.  Prony  derived  residues  and  poles  provide 
an  abbreviated  storage  mode.  However,  no  physical  significance  is  attached  to  the  Prony 
derived  coefficients  and  arguments.  They  are  used  only  for  reconstruction  of  the  original 
data.  As  an  example,  a 180-point  representation  of  a pattern  could  be  represented  by  6 
to  10  data  points. 

The  data  compression  attributes  of  Prony’s  method  can  also  be  used  to  develop 
simpler  equivalent  source  distributions  for  complex  antennas.  Currently,  this  equivalent 
sources  concept  is  being  pursued  within  the  method  of  moments  and  geometric  theory 
of  diffraction  efforts  at  NOSC. 
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A.l  INTRODUCTION 


Angle  Domain  Pattern  Analysis  Program  (ADPAP)  performs  either  a Prony  or  a 
Fourier  series  analysis  on  radiation  patterns  for  data  supplied  in  equal  degree  increments 
for  any  7r  interval  in  the  range  0 < 0 < 27r.  Both  magnitude  and  phase  of  the  radiation 
pattern  are  required  for  correct  results.  The  Fourier  series  expansion  may  be  applied  to 
either  the  real  or  the  imaginary  part  of  the  pattern. 

The  shifting  theorem  may  be  applied  to  both  the  Fourier  and  the  Prony  analysis 
to  move  the  apparent  phase  center  of  the  radiation  pattern.  One  shift  in  both  X and  V 
directions  may  be  applied  before  selection  of  the  n interval  for  processing.  Multiple 
shifts  along  the  assumed  array  axis  may  be  applied  after  selection  of  the  v interval  and 
prior  to  the  initiation  of  the  Fourier  or  Prony  analysis.  The  shift  can  be  specified  in 
feet  or  in  metres. 


A. 2 DATA  INPUT  FORMAT 

ADPAP  uses  a free  format  input  all  data  items  on  a single  line  are  separated 
by  commas. 

Free  format  input  is  accomplished  by  suspending  the  standard  storage  control 
table  element  which  normally  checks  for  record  type  when  reading  from  a file.  For 
the  NOSC  UNI  VAC  1110  machine,  this  is  accomplished  at  time  of  compilation  by 
including  the  following  commands  in  the  run  stream  between  the  compilation  and  map: 

("COPY  ,S  N*FTNLIB. 

("ASM  ,L1  F2SCT  .F2SCT 
FSSCT33  .111  .0 
END 


A. 3 DATA  CARDS 

ADPAP  uses  a data  input  format  similar  to  NEC  Method  of  Moments' . Each 
line  of  data  begins  with  a two-letter  mnemonic  which  determines  how  the  rest  of  the 
data  on  that  line  will  be  interpreted  bv  the  code.  Each  data  set  must  begin  with  at  least 
one  comment  card.  Any  number  of  comment  cards  may  be  used,  but  the  last  one  must 
be  a CE  comment  card.  (If  only  one  line  of  comments  is  required,  then  that  card  must 
be  a CE  card.)  The  order  of  the  other  cards  is  unimportant  except  that  the  XQ  causes 
execution  of  the  program  for  the  parameters  set  on  all  previous  cards  back  to  the  last 
XQ  card,  if  there  is  one.  Use  of  the  NX  card  re-initializes  all  program  parameters, 
allowing  introduction  of  a new  problem  within  the  same  data  stream. 

All  data  cards  except  for  the  comment  cards  have  a two  letter  mnemonic  followed 
by  two  integer  fields  and  two  floating  fields.  One  other  special  case  also  differs  from 
this  format  which  is  explained  under  the  Pattern  Data  Command  (PD)  card. 


1 . Naval  Ocean  Systems  Center  Technical  Document  1 16.  Numerical  Electromagnetic  Code  (NEC) 
Method  of  Moments,  by  GJ  Burke  and  AJ  Poggio,  18  July  1977. 


COMMENT  CARDS  (CM,  CE) 


CM  - COMMENTS  - 


CE  - COMMENTS 

NOTES 

Any  number  of  CM  cards  may  be  used. 

The  last  comment  must  be  a CE  card. 

For  only  one  line  of  comments,  use  a single  CE  card. 
A maximum  of  one  CE  card  may  be  used. 

ALL  comment  cards  must  appear  first  in  the  data  set. 
At  least  one  comment  card  must  follow  an  NX  card. 

END  of  RUN  (EN) 


PURPOSE:  To  indicate  to  the  program  the  end  of  all  execution. 

CARD:  EN 

FOURIER  ANALYSIS  (FA) 

PURPOSE:  To  initiate  a Fourier  analysis 

CARD:  FA,  NTERMS,  NSAMPS,  START 

PARAMETERS: 

Integers 

NTERMS  Desired  number  of  Fourier  coefficients 

NSAMPS  — Number  of  numerical  integration  increments  used  for  Euler 
Equations 

Decimal  Numbers 

START  Initial  angle  of  180°  of  patterns  used  in  analysis  in  degrees. 


NOTES: 


The  range  of  source  locations  is  (NTERMS)  meters. 

The  possible  source  locations  are  at  S = ± n(X/2),  n = 0,  1,2, NTHRMS. 

NTERMS  greater  than  0;  REAL  (E(0))  is  processed 
NTERMS  less  than  0;  1MAG  (E(0))  is  processed 

Only  180°  of  pattern  may  be  analyzed  at  a time.  An  180°  sector  of  the  pattern 
may  be  used  beginning  with  the  angle  START  (in  degrees).  Symmetry  of 
pattern  is  always  assumed. 

NEW  PROBLEM  (NX) 

PURPOSE:  To  signal  the  end  of  one  problem  and  the  beginning  of  an  entirely  new  problem 

CARD:  NX 

NOTES: 

This  card  will  re-initialize  all  control  parameters  and  arrays. 

The  card  that  directly  follows  the  NX  card  must  be  a comment  card. 

PRONY  ANALYSIS  (PA) 

PURPOSE:  To  initiate  a Prony  analysis 

CARD:  PA.  NPOLES.O.  START 

PARAMETERS: 

Integers 

NPOLES  Desired  number  of  poles 
0 Blank 

Decimal  Numbers 

START  Initial  angle  of  180°  of  pattern  used  in  analysis  in  degrees 

NOTES: 

The  number  of  pattern  data  samples  is  equal  to  twice  the  number  of  poles 

Only  180°  of  pattern  may  be  analyzed  at  a time.  Any  180°  sector  of  the 
pattern  may  be  used  beginning  with  the  angle  START  (in  degrees).  Symmetry 
of  pattern  is  always  assumed. 

PATTERN  DATA  (PD) 

PURPOSE:  To  specify  the  sampling  parameters  and  data  of  the  radiation  pattern 

A4 


w 


CARD: 


(Format  A) 

PD,  NF,  ICVR,  F,  ELW 
(Format  B) 

0(1)  E(l) 


Phase  ( 1 ) 


0 (NF)  E (NF)  Phase  (NF) 

PARAMETERS: 

Integers 

NF  - Number  of  field  points 
ICVR  - Special  processing 

IF  = -1 , Solves  for  amplitude  of  an  elementary  dipole 
IF  = 0,  point  sources 

IF  = +1,  Solves  for  amplitude  of  X/2  dipoles 

Decimal  Numbers 

F - Frequency  in  Hz 

ELW  - If  ICVR  = -1,  effective  length  in  wavelengths  of  elementary  dipole 
(I)  - Pattern  angle  of  sample  I 
E (I)  - Magnitude  of  field  of  sample  I 
Phase  (I)  Phase  of  field  of  sample  I 

NOTES: 

The  first  card  follows  format  A with  the  definitions  given  above.  The  next  NF 
cards  have  format  B giving  the  azimuth  pattern  angle,  field  strength  magnitude 
and  phase  as  indicated. 

PRINT  CONTROL  (PT) 

PURPOSE:  To  control  printing  of  pattern  data 

CARD:  PT,  IPT1.IPT2 

PARAMETERS: 

Integers 

IPT1  - = 0,  no  output  of  original  pattern  data 

= 1,  prints  original  pattern  data  in  equal  cosine  steps 

IPT2  - = 0,  no  output  of  generated  pattern 

= N,  N greater  than  0,  computes  and  prints  pattern  using  Fourier 
coefficients  for  0 = 0,  180°  in  N equal  steps 
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* . 


Decimal  Parameters 
No  parameters 

NOTES: 

To  erase  a previous  PT  specification,  either  use  a new  PT  spec  or  use  PT.  0, 
or  use  NX  to  re-initialize  all  arrays. 

If  no  PT  is  specified  in  a given  run,  no  pattern  data  will  be  output. 
SHIFTING  THEOREM  (SF) 

PURPOSE:  To  multiply  the  pattern  data  by  a phase  shift 


CARD:  SF,  ITER,  UNITS.  X,  Y 

PARAMETERS: 

Integers 

ITER  less  than  0:  resets  pattern  to  original  pattern 

= 0;  shifts  phase  center  to  location  (X,  Y)  greater  than  0;  number 
of  iterations  of  shifting  theorem  applied  along  array  axis 

Decimal  Numbers 

X ITER  = 0:  X = shift  along  X axis 

ITER  greater  than  0:  X = initial  shift  along  array  axis 

Y ITER  = 0;  Y = shift  along  Y axis 

ITER  greater  than  0.  Y = increment  of  shift  along  array  axis 

NOTES: 

For  ITER  greater  than  0,  the  factor  exp  (ikXcos  9)  is  applied  to  the  pattern 
after  incrementation  in  equal  cosine  steps  and  after  rotation  of  the  pattern  by 
START  (specified  by  the  FA  and  PA  commands). 

For  ITER  = 0,  the  factor  exp  (-ik(Xcos  + Ysin))  is  applied  to  the  initial  pattern. 
For  ITER  less  than  0,  all  previous  SF  specs  are  erased  and  the  original  pattern 
will  be  processed. 

The  array  axis  is  the  positive  X axis  or  the  line  defined  by  the  angle  START  on 
the  FA  and  PA  cards. 

TWO  SF  cards  may  be  used  consecutively  before  the  XQ  so  that  both  shifting 
options  can  be  used. 

EXECUTE (XQ) 

PURPOSE:  To  cause  program  execution  at  desired  points  in  the  data  stream. 


CARD:  XQ 

NOTES: 

The  XQ  card  is  required  to  cause  execution  of  the  program. 
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4 
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12 
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15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 
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28 

So9 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

n 

50 

51 

52 

53 

54 

55 

56 


C 4 4 4 4 

c * * * « 
C 4 * * 4 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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ANGLE  DOMAIN  PATTERN  ANAIVSIS  PROGRAM* 


fOURIER  SERIES  EXPANSION 

DEVELOPED  BY  J.  ROCKWAY  AND 
APRIL  1978 

ANGLE  DOMAIN  PRONV  


C.  LOGAN 


C 

500 


501 


502 

100 

1 

101 


ADAPTED  EROM  L L L CODE  BY  J.  C.  L06AN 
APRIL  1978 

COMPLEX  E <1001 ,EP<360  It  A { 100) ,ALPMA<  100) ,f PP  < 360> ,EO(  100) .2 , SUM 
DIMENSION  CI(51).CJ<50),S<51).AMP(100),PH(100) 

DIMENSION  AC<10C),AD(100),AP(360>,SCAIE(2>,EI(100>.Ey( 100) 
CHARACTER  BLANK/'  ' / , SL ASH/'/ '/ , C ARD (78) 

CHARACTERS  IR<?6  >,  I T (10) 

DATA  IT/2HCM,2MCE,2MPD,2HP»,2M»8,?HNK,2HEN,2HfA,2MPT,?MSf/ 

DATA  Pl/3. 14  15927/ 

DATA  SCALE/1 .,0.304801/ 

DEEINE  EILE  9 ( S D E ) 

INITIALIZE  OATA 
I PR  ONY  * 0 
I EO  UR*C 
IRE  AL*1 
START*0, 

SX=0.0 
SY*0.0 
S A*  0 . 

DSA  *0. 

ELS. 

ELH-O. 

1 C U R = 0 
I PT 1 =0 
1 PT  2 *0 
ITE  R=0 

DO  501  1*1,360 

A P ( I)*0. 

E P ( I >*CMPLX (0. .0 . ) 

EPP  (I)*EP(I > 

CONTINUE 
DO  502  1*1,100 
E (I  >*CMPLX(0.,0.  ) 

AlPHAU  ) = CMPL  X(0  . ,0.  ) 

A ( I >*ALPHA( I ) 

AD( I )*0. 

AC ( I )=C. 

CONTINUE 

WR1TE(6,100) 

EORMAT ( 1 HI , TIB. ' • 44  44  angle  DOMAIN  PATIFRN  ANALYSIS  PROGRAM  ' 

1'44444  '///TY6,C0"MENTS'/T36,' '//) 

RE  A 0 ( 5 , 1 01  ) 1R 
E OR  MAT  ( 36  A2  ) 

VR1TE(6.102>  ( 1R( I) , I *2 ,36  ) 


b-: 


57 

102 

38 

59 

60 

61 

103 

6? 

63 

200 

64 

65 

66 

67 

68 

104 

69 

70 

5 

71 

72 

105 

73 

74 

75 

6 

76 

77 

78 

7 

79 

8 

80 

i 

81 

82 

83 

84 

85 

C 

86 

£ i 

87 

c 

88 

C CM 

89 

C CE 

90 

C PD 

91 

C PA 

92 

C X Q 

93 

C NX 

94 

C EN 

95 

C FA 

96 

C PT 

97 

C SF 

98 

C 

99 

100 

106 

101 

103 

104 

’Q5 

106 

107 

108 

109 

110 

115 

111 

ill 

C 

FORMATC5X.35A2) 
i f < r r ( i ) .eq.itci ) > so  to  i 

IF(IR<1).E«.IT<2))  60  TO  200 
WRITEC6.103) 

FORMAT<2X,'READ  ERROR  IN  INPUT  DATA  - STOP!') 
STOP 

ITEMP1=0 
ITE  MP2  = 0 
TEMP1=0. 

TENP2-0* 

READC5.104)  1R(1),CARD 
FORMATE A2 *78 Al) 

C AR  D ( 1 ) =6  LANK 

I F < CARDC78) .EO. blank ) 60  TO  6 
WRITEC9.105)  CARD 
FORMATC78A1  ) 

REAOC5.105)  CARD 
GO  TO  5 
DO  7 1=2,78 

IF (CARD(l) ,E0. BLANK ) CAR0(1)=SLASH 
IF < CARD ( 1 ) .EQ  . SLASH)  GO  TO  8 
CONTINUE 
CONTINUE 

WRITE(9,105)  CARD 
END  FILE  9 
REWIND  9 

REAO(9,«,ERR*1000)  ITEMP1,ITEMP2,TEMP1,TEMP2 
REWIND  9 

COMMANDS  : 

COMMENTS  1 
COMMENTS  2 
PATTERN  DATA  3 
PRONY  ANALYSIS  4 
EXECUTE  5 
NEW  PROBLEM  6 
END  EXECUTION  7 
FOURIER  ANALYSIS  8 
PRINT  PATTERNS  9 
SHIFTING  THEOREM  10 

MR I TEC 6, 106)  IRC1 ), IT EMP1.IT EMP2, TEMPI  ,TEMP2 
FORNAT<//2x,3H***.A2,214,2E15.5) 

I F < IRC1) .EO.ITCI ) ) 60  TO  30 
IF  < I R ( 1 > .EQ.ITU)  ) 60  TO  4 0 
IF(IR(1).E9,IT(5>)  GO  TO  50 
I F C IRC1) .E0.ITC8) ) 60  TO  800 
IF< IR(1) .E8.ITC9)  ) GO  TO  900 
IFC IRC1) .ES.ITCD)  ) GO  TO  5 00 
1 F < 1RC1) .E6.1TC10))  60  TO  600 
IF<IRC1).E«.ITC7>)  GO  TO  700 
WRITEC6.115) 

F0RMATC///2X, 'COMMAND  IS  UNDEFINED  - STOP!') 
STOP 


READ  IN  PATTERN  DATA 


114 

C 

115 

30 

116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 
127 

31 

128 

129 

32 

130 

131 

132 

133 

134 

33 

135 

136 

34 

137 

138 

139 

140 

141 

35 

142 

143 

C 

144 

C 

145 

C 

146 

40 

147 

148 

149 

150 

151 

152 

C 

153 

c 

m 

c 

600 

156 

610 

157 

158 

159 

160 

620 

161 

162 

163 

630 

164 

165 

166 

W 

169 

C 

170 

C 

NF=ITEMP1 
ICUR=1TEMP2 
F=TEMP1 
ELN=TEMP2 

1 r < F ,£Q.  0.)  F=299.8E*06 
W = 2 .998  E ♦ 08  / F 
PK=6.2831853/W 
00  31  1 = 1 , NF 
REAMS,*)  AP(1),X,V 
YRAD=0.0174533*Y 
XR=X»C0S ( TRAD) 

YI*X*SIN(YRAD) 

E0«  I)  = CPPLX(XR,  Y1  ) 

IFIITEMP2J  34,200,32 
IFtELU  ,£Q.  0.)  E LW  =0  ,5 
EL=ELW*W 
DO  33  1=1, NE 
APR=AP(I)*0. 0174533 

FTR=-SIN(APR)/<COS<0.5»PK«EL*COS<APR )>-C0S(0.5*PK«EL)) 
£Ofl>=EO(I)«CMPLX(0.,FTR/60.) 

60  TO  200 

IF(ELU  .EQ.  0.)  ELR =0,125 
EL=ELW«W 
00  35  1=1 ,NF  , 

APR=AP<I )*0. 0174533 
FTR=-30.«PX*EL*SIN(APR) 

EO< I)=EO(I)*CMPLX (0,, 1/ETR ) 

GO  TO  200 

SET  PARAMETERS  FOR  PRONY  ANALYSIS 

NP0LES*ITEMP1 
START=T£HP1 
NSAMPS=2*NP0LES 
IPR0NY=1 
1 FOUR=0 
GO  TO  200 

SET  UP  PARAMETERS  FOR  SHIFTING  THEOREM 

IF  < I TEMP  1 ) 630,520,610 
ITE  R = ITEMP1 

DSA=TEMP2*SCALE(ITEMP2+1 ) 

SA* TEMPI* SCALE! I TEMP2*1> 

GO  TO  200 

SX= TEMPI* SC  ALE O T EMP2 *1 ) 

SY=TEMP2*SCALE(ITEMP2+1) 

GO  TO  200 

SA  = 0 • 

ITE  R = 0 
OSA  =0. 

SX=0.0 
SY=0.0 
60  TO  200 

- SET  UP  PARAMETERS  FOR  FOURIER  SERIES 


I f OUR* 1 
1 PR  ONT  *0 

NTERMS*aBS(ITeNP1  > 

NSAMp$*ITEMP2 

START*TEMP1 

IFUTEMP1  .LT.  0)  IRE  AL  r 0 
GO  TO  200 

PRINT  CONTROL 

IPT 1=ITEMP1 
1PT2*ITEMP2 
60  TO  200 

PROGRAM  EXECUTION 

CONTINUE 

IMS*  .ED.  0.0  .AND.  ST  .EQ.  0.0)  GO  TO  693 

APPLY  SHIfTING  THEOREM 

DO  691  1*1. NF 

APR*AP<I )«0. 017*533 

AR6*-PK«  < SX«COS<  APR )*ST«S1 N(APR ) ) 

2*CMPLX<0..ARG) 

EPP ( I>*EO( J ) *CEXP (2 ) 

CONTINUE 
GO  TO  692 
CONTINUE 
DO  69*  1*1, NF 
EPP  <I)*EO(I  ) 

CONTINUE 

ITER  * R OF  ITERATIONS  OF  SHIfTING  T HE  OR M 
DO  999  IST£P*1,ITER 
WRITE  <6.127) 

F OR MATI///T36,3(' *••••')//) 

WR1TE<6.126)  SX.ST 

FORMAT<///5X,THE  phase  center  is  LOCATED  AT  (', 
1M0. *,',',E10. *,'»', 3 X, 'METERS'//) 

WRITE<6,117)  F.u 


format<//5x.'freouency  *'.E15.*,  , 

1'  H2'//5X, 'WAVELENGTH  * »F  1 5.  * , 

-select  ?pi  range  to  be  processed 

IFtSTART  .EQ.  0.)  GO  TO  10 
FIN*START*HO. 

I F < F I N .6T.  360.)  FIN *F IN- 360. 

DO  11  1*1, NF 

I F ( START  .EQ.  AP<D)  JS*I-1 
I F < F IN  .EQ.  AP ( I ) ) J F *1  - 1 
CONTINUE 
ISK  IP*JF-JS 
IF(ISKIP)  12,13,1* 

CONTINUE 

L1M*NF-J$«1 


METERS'//) 


228 

H 3 

15 

231 

232 

233 

1 6 

234 

235 

14 

236 

237 

238 

1 7 

239 

240 

13 

241 

116 

itf 

10 

244 

245 

18 

1 9 

248 

ill 

C 

C 

251 

c 

252 

253 

254 

255 

256 

257 

60 

258 

c 

259 

c 

260 

c 

i*3 

c 

263 

264 

265 

266 

267 

268 

70 

269 

80 

270 

271 

272 

90 

273 

274 

C 

275 

c — — 

276 

c 

277 

278 

279 

280 

640 

281 

282 

125 

283 

284 

690 

00  15  1*1  .LIP 
EP<I)*EPPf JS*I> 

CONTINUE 
00  16  I * 1 , J F *1 
E P ( I*LIP)=EPP(I> 

CONTINUE 
60  TO  19 
CONTINUE 
00  1 7 1*1 ,NF 
EP<  I)*EPP(1*JS) 

CONTINUE 
60  TO  19 
WRJ  TE(6,116) 

F0RPAT<2x, 'error  IN  START  SPECIFICATION-STOP*') 
STOP 

CONTINUE 
00  18  1*1, NF 
EP(I)*EPP<I) 

CONTINUE 

CONTINUE 

0ELC0S=-2 ./ (NSAPPS-1 ) 


SET  UP  AN6LES  FOR  EQUAL  COS  SPACING 


AC ( 1 )*1  . 

AD( 1 )*57 ,2957795 *ACOS (AC (1 )) 

00  60  I=2,NSAPPS 

AC  < I)=AC  < 1-1 )*DELCOS 

ADI  I ) — 5 7 .2957795* ACOS (AC (I  )) 

CONTINUE 


DETERMINE  PATTERN  IN  EQUAL  COS 
0 V LINEAR  INTERPOLATION 

E (1 > = EP<  1> 

E ENSANPS )*EP(NF ) 

00  90  J=2.NSA«PS 
DO  70  1*2, NF 

JN*  I 

1F(AP(1> .GT.ADIJ ) ) 60  TO  80 

CONTINUE 

CONTINUE 

FCTR*(AP(JP)-AD(J  ))/(AP(JP)-AP(JP-1)  ) 

£(J)*EP(JP)-<EP<JP>-EP(JP-1))*FCTR 

CONTINUE 

IF(  ITER.EQ.O)  60  TO  6 90 

APPLE  SHIFTING  THEORE*  ALONG  ARRAY  AXIS 


00  840  I=1,NSAPPS 
Z*CPPLX(0.,PK*AC ( I ) • S A) 

E(1)*E(I)*CEXP(Z) 

CONTINUE 
WRITtC6,125>  SA 

E0R"AT(///5X,  THE  PHASE  CENTER  HAS  BEEN  SHIFTED  BY' 
1 F12  ,3X  , 'PETERS  ALONG  ARRAY  AXIS'//) 

CONTINUE 


+>  *• lMW*l*V*OQ KMiMXA SMXjtTV f\jr\>fN>f\/l\) rv/rofM tv -*-*-*-*-* •*• -*  -* OOO  OO OO OO 0*0<0  'OO'© «« 'O'©  <OOMl QNMI 

-^0>OQB'>40^#'04f\i-*0'0  00^40v^#*L^>j  -+0'0(»’^<* V***Ui*J-«0<OaD'40t/' ^UR'V-AQOOD N>^W'V-*0'OOI ^9>UI 


c 

c 

c 


91 

107 


108 

901 

109 

C 

c 

c 

c 

c 


95 

110 

111 

112 


115 

801 

C 

C 

C 

c 

c 


806 

807 

808 
809 


IM  I P T 1 .EQ.  0)  (0  TO  901 

PRINT  PATTERN  IN  EQUAL  COS  (THETA) 

DO  91  J * 1 .NSAMPS 
APP ( J)*C ABS(E( J) > 

PM(J)*180.*9TAN2<A1NAG(E<J)),REAL(E<.I>>>/PI 

CONTINUE 

NRITE(6,107) 

FORMAT (//2X. 'PATTERN  DATA'//7X, 'THETA', 2X, 

1 'C OS(THETA), 10X, 'R EA L.IOX.'IMAG.I OX, ‘mAG'.IOX, 'PHASE'/) 

WR1 TE(6,106)  (ADd),AC(I  ),E( D.ARPd  >,PH(I), 1*1 .NSAMPS) 
FORmAT(5X,F7.3,F12.5,3£1S.5.F12.3> 

CONTINUE 

HRITE(6,109)  DELCOS 
F0RMAT(//5X,'DELC0$  *',E15.5) 

IFCIPRONY  .EQ.  0)  60  TO  801 

8E6IN  PBONY  PROCESS 

CALL  PR0NT(0ELC0S ,E ,A  ,AL  PH  A , NP0  L E S ,N SAMPS) 

CORRECT  RESIDUES  FOR  SERIES  NOT  STARTING  AT  EXP(0) 

NORMALIZE  POLES  TO  UNITS  OF  DISTANCE 
DO  95  1*1 »NP0L  E S 
A (I  )*A(I )*CEXP(-ALPHA (1 ) ) 

Z*A (1) 

PH(I)*180«*BTANZ(AIMAG(Z).REAL(Z))/PI 
ALPHA(1)*ALPHA(I ) /PK 
AMP ( I )*C A0S ( A(  I ) ) 

CONTINUE 

WRITE(6,110) 

F0RMAT(//2X, -COMPL EX  PRONT '/ / ) 

NRITE(6,111)  START, NSAMPS.NPOLES 

FORMAT(//5x, 'START  * '.F8.3.'  D E 6RE E S ' / / 5X , 'NO . OF  SAMPLES  *', 
TI5//5X, 'NO.  OF  POLES  * .15//) 

I F ( IC UR  .EQ.  -1)  WRITECo.l 30)  ELM. EL 
I F ( ICUR  .EQ.  1)  WRITE  (6,12  9)  ELW.EL 
WRITE(6,112>  , ... 

F OR MAT(//16X. 'RESIDUES, 50 X,'P0LES/12X, 
1'REAL',12X,'IMAG',12X,'MAG',10X,'PHASE',10X,'REAL,12X,'IMAG'/) 
WRI TE (6.113)  (1 ,A (I ) , AMP (I ),PH( I ) ,AL  PHA( I) ,1*1 ,N POLES) 
F0RMAT(2x,I2,3E15.5,F15.3.2E15.5) 

I F ( I FOUR  .EQ.  0)  GO  TO  990 

BE6IN  FOURIER  SERIES  ANALYSIS 

IF  IREAL*0, EXPAND  IMA6(E>  IN  FOURIER  SERIES 
IF  IREALd, EXPAND  R E A L ( £ ) IN  FOURIER  SERIES 
IFdREAL  .EQ.  0)  GO  TO  807 
DO  806  I *1  .NSAMPS 
EX(I)*REAL(E(I)) 

ET(1  >*AIMAG(E( I) ) 

GO  TO  809 

DO  808  1*1. NSAMPS 

EY(I)*REAL(E(I>) 

EX(1)*AIMAG(E(I>) 

CONTINUE 


370 

371 

372 

373 

374 

375 

376 

377 

378 

379 

380 

381 

382 

383 

384 

385 

386 

III 

389 

390 

391 

392 

393 

394 

395 

396 

397 

398 


CALL  FOUR(-DELCOS,AC,EX,CI.CJ,NTERNS.NSAHPS) 

C DETERMINE  SOURCE  STREN6TH  AND  PHASE 

If ( 1REAL  • E Q • 0)  GO  TO  810 
CALL  fOURO(-OELCOS,ET  ,C0,NSANPS) 

PH(  1 ) = 8TAN2(C0,C I (1 ) ) 

A«P(1)*CI(1)/COS(PH(1 )) 

PH(  1 ) =57 . 29  57795  *PH ( 1 ) 

DO  805  1*1,NTERNS 
PH(I*1)=8TAN2(-CJ(I),CI(I+1)) 

ANp(I«1)=CI<I*1)/COS(PH(I*1)) 

PH(  1*1  )=57.2957795*PH (1*1 ) 

NRITE(6,118) 

FORNAT(//2x,' FOURIER  COEFFICIENTS  FOR  REAL(E)'//) 

GO  TO  812 

CALL  FOUR0(-DELCOS,ET  .CO.NSAhPS > 

PH(1)=BTAN2(C 1(1), CO) 

A N p (1)  =C  I ( 1 ) /SIN  ( PH  ( 1 )> 

PH(1)=57,2957795»PH(1 ) 

DO  dll  I=1,NTER*S 
PH( l4l)=BTAN2(CI  ( 1+1  ) ,CJ (I  )) 

ANP<l«1)=CI(I*1)/SIN(PH(l*1>) 

PH<I*1)*57.295  7795»PH(M1) 

WRlTE(6,126) 

fornat<//2x,' — Fourier  coefficients  for  inag(E)'//) 

CONTINUE 

DETERMINE  SOURCE  LOCATIONS 
DO  804  I*1,NTERHS*1 
S(I)=0.5«(1-1)*W-SA 
WRITE(6,119)  start, nsanps.nterns 
FORNAT(//5X, 'START  = '.Fg.J  ' 0 E GRE E S ' / / 5* , 

1 ' NO  • OF  SANPLES  = '.I5//5X,  NO.  OF  TERNS  * ',I5//> 

1 f ( I C UR  • E 9 • -1)  WRITE(6.130)  E L W , E L 
FORNAT( //5X , 'CURRENTS  OF  ELENENTAAT  RADIATORS  : 

1'TH  =',f 12,4,  LANBDA  (OR  .F12.4,'  NETERS)  > 

IF ( 1 C UR  ,EO.  1)  WRITE (6,129)  ELK. EL 
FORNAT(//5X, 'DIPOLE  CURRENTS  : TOTAL  LENGTH  = 

1'  LANBDA  (0R',F12.4,'  NETERS)') 

WRITE(6,120> 

F0RNAT(//!x,'N',4X,'S(NETERS),6X, 

I'l (N) ',12X,'J (N) ' ,20X ,'ANP (N)',8x .'PHASE '/) 

WRI TE (6.121)  S(1),CI(1),ANP(1),PH(1) 
FORNAT(3x,'0',F12.4,1E15,5,25x,E15.5,F12.3) 

WRITe(6.122)  (I ,S  (I *1  ).C  I ( 1*1 ) . C J (1)  ,ANP ( I ♦ 1 ) ,PH ( 1*1 ) ,1*1 ,N TERNS ) 
FORNAT(2X,I2,F12.4,2e15,5,10x,1E15.5,F12,3) 

I F ( IPT2  ,EQ.  0)  60  TO  990 

CALCULATE  AND  WRITF  PATTERN  FRO«  FOURIER  COEFFICIENTS 
WRI TE(6,123) 

FORNAT(//2X, 'PATTERN  FRON  FOURIER  SE R I E S '/ / 7 X , ' T HE T A ' , 1 OX , ' R E AL ' , 
110X,'INAG,10X,'NA6', 10*. 'PHASE'/) 

XB=18Q./IPT2 
DO  802  I=1,IPT2*1 
BE  T A * ( I — 1 ) • X§ 

BETAR=BETA«PI/180. 

SUN*CNPLX(0.,0.) 

DO  803  J=1,NTERNS»1 

AR6=PK*S(J)*COS(BETAR)*ph(J)»0.0174533 
803  SUN  =SUN«ANP(J ) *CNPL  X ( COS  ( A RG ) ,S I N ( AR G ) > 


805 

118 

810 


811 

126 
812 
C 

804 

119 

130 

129 

120 

121 

122 

C 

123 


EFFECTIVE  LENG' 


i F 1 2 ,4 , 


1F(8ETAR  . E 0 • 0.  .^R.  BETAS  .E8.  180.)  60  TO  711 

likW  •<COS?0.?*PK*EL»COS <BETAR))-C0S<0.5«pK»EL))/SINC8ETAR> 
SUN=SUH*CMPLX(0., ARG) 

60  TO  730 
SUH=CNPLX<0.,0.) 

SUH  A *0  . 

PHA  S£ *0 • 

60  TO  740 

AR6*30.*PK*EL*S1N(BETAR) 

SUH=SUN*CHPtX (0. , ARG) 

CONTINUE 

SUMA-CABS(SUM) 

PHASE *180.»8TAN2< AIHA6<SU8 ),REAL (SUN )>/PI 
HRITE(6,124)  BETA,  SUM  .SUHA.PHASE 
FORHAT(5x,F7.3,3E15.5 ,F12.3) 

CONTINUE 
SA=SA*DSA 
CONTINUE 
GO  TO  200 
«RITE(6,114) 

F0RNATC2x,'REA0  ERROR  ON  UNIT  9 - STOP!  ) 

STOP 

END 


i 

J 

4 

5 

6 

7 

8 
9 


B-9 


FUNCTION  BTAN2<X1,I2) 
BT*n2*0.  WHEN  xl  = x 2 -o . 
IF<>r1>  3,1,3 
1FCX2)  3.2,3 
B TA  N2r0 • 

RETURN 

8TAN2  = AT  AN2  C XI ,X2  ) 

RETURN 

ENO 


subroutine  crout  <C,  X,  »,  NCOLSC.  N COL  ST) 

FORTRAN  CROUT 

COUPLE*  SUM.H I 6H.S UM 1 ,0 ET£RM 
CONPLEX  A,  C,  X,  T 

OINENSION  A ( 50 . * 1 > , C (50,50),  *(50.1  ),  *(50.1  ).  INDE*(50> 
CROUT  REDUCTION 


THIS  IS  A CROUT  REDUCTION  PROGRAH  WHICH  CAN  EITHER  SOLVE  FOR  TmE 
SOLUTION  OF  A NATRIX  EQUATION,  OR  AS  A SPECIAL  CASE,  IT  CAN  CONPU 
INVERSE  OF  A SPECIFIED  N BY  N "ATRIX 

THE  TECHNIQUE  IS  CLEARLY  EXPLAINED  ON  PAGES  4 ? 9 TO  435  OF  THE  BOOX 
F.P.  HILDEBRAND.  INTRODUCTION  TO  NUMERICAL  ANALYSIS,  PUBLISHED  BY 
MCGRAW  HILL  IN  1956 

THIS  PROGRAM  SOLVES  THE  MATRIX  EQUATION  HX  = B.  THE  H AND  THE  B 
MATRICES  ARE  SPECIFIED  PARAMETERS,  AND  THE  PROGRA"  CALCULATES  TH 
MA  TR I X . 

H IS  an  N PY  N MATRIX  OF  NUMBERS,  B IS  AN  N BY  M MATRIX  OF  NUMB 
AND  X IS  AN  N BT  M MATRIX  OF  NUMBERS  TO  BE  DETERMINED.  IF  THE  IN 
OF  THE  H MATRIX  IS  DESIRED,  ONE  SETS  THE  B MATRIX  EQUAL  TO  AN  N 
IDENTITY  MATRIX,  THUS  X BECOMES  THE  INVERSE  MATRIX. 

THE  A MATRIX  SPECIFIED  IN  THE  ARGUMENT  OF  THE  SUBROUTINE  IS  AN 
AUGMENTED  MATRIX  CONTAINING  THE  N BY  N MATRIX  H IN  THE  FIRST  N BY 
LOCATIONS,  AND  IS  AUGEMENTED  IN  THE  N ♦ 1 TO  N ♦ M COLUMNS  BY  THE 
MATRI X B. 

NN  IS  BY  DEFINITION  EQUAL  TO  N ♦ M. 

THE  SOLUTION  MATRIX,  X,  IS  STORED  IN  THE  FIRST  N BY  M (ROWS  p Y COL 

ELEMENTS  OF  THE  AUGMENTED  MATRIX  A,  THE  ORIGINAL  H AND  B MATRICES 
DESTROYED,  AND  SHOULD  BE  SAVEO  ELSEWHERE  IF  THEY  ARE  TO  BE  USED  AF 
USING  THI S ROUTINE . 


COMMENT SOLVES  C X = Y FOR  * C=COEFF 

COMMENT COPIES  C AND  Y INTO  A MATRIX 

N = N C OL  SC 
M=NC0LSY 
NN=  N*M 

DO  2 QQ  I Q 1 = 1 , N 
DO  £00  I Q 2 = 1 , N 

200  AUQ1.I02)  = CCIQ1.IG2) 

DO  210  I Q 3 = 1 ,N 
DO  210  104=1, m 

210  A(IQ3,IQ4*N)  * Y(!Q5,IQ4) 

ABJC  - 1 .E -20 
J 1 =N-1 
JA  =N*1 


C=COEFF  MATRIX 


tO  30  1 1 1)  N 
INDEX(I)=I 


oo  roo  j=i,nn 

0 0 8 00  II=1.N 

sum=o.o 

1=INDEX(II) 

IF(I1-J)33,34,34 

THE  SECTION  AFTER  STATEMENT  33  IS  EQUIVALENT  TO  FORMULA  10.4.5  AN 
10.4.6  OF  N1L0EBRAN0 

33  If  <11  - 1)  9000,9200,9000 
9000  LLLL  =11-1 

DO  9100  X = 1 , L L LL 
IPPP  = INDEXOO 

9100  SUM  = SUM  ♦ A( I ,X >» A( IPPP, J > 

9200  A ( 1 , J ) = (A(ItJ)-SUM)/A(I,II) 

GO  TO  800 

THE  SECTION  AFTER  STATEMENT  34  IS  EQUIVALENT  TO  FORMULA  10.4.4 
OF  HILDEBRAND 

34  IF  (j  - 1)  8000,8200,8000 
8000  LLLL  = J - 1 

DO  8100  X=1,LLLL 
IPPP*INDEX<K) 

8100  SUM=SUM*A{I,K)*A( IPPP , J ) 

8200  All  ,J)  = A(I,J>— SUM 
800  CONTINUE 


IF<J-NU1,700.700 
41  L = I NDE  X ( J ) 

THIS  SECTION  SHIFTS  AND  REORDERS  THE  COLUMNS  AND  ROWS  ( THEY 
RESHUFFLED  AT  THE  END  Of  THE  PROCESS  TO  BE  PUT  IN  THE  ORIGINAL  0 


X A = L 

HIGH=A(L«J) 

X 7=  0 

00  35  I = J , J Z 
JC=I*1 
L*INDEX< JC) 

IF  (CABS(HIGH)  - CABS(A(L,J>)>  36,35,35 
36  NIG  H*A ( L , J ) 


1U 

m 

117 

118 

119 

120 
121 
122 
123 
12* 
125 

w 

128 

129 

130 

Hi 

13! 

13* 

135 

136 

137 

138 

139 
1*0 
1*1 
1*2 
1*3 
1** 
145 
1*6 
147 
1*8 
1*9 

150 

151 

152 

153 
15* 

155 

156 

157 

158 

159 

160 
161 
162 
163 
16* 

165 

166 


W 

170 


KA*L 
KZ  = 1 

35  CONTINUE 

IF  (CABS (HI6H)  - ABJ C)  31,31,3200 
31  WR1TE<3,32)  ABIC 

32  FORMAT  ( 2X  , 2*HP I VOT  ELEMENT  LESS  THAN  ,E20.8) 
3200  00  37  K = 1 ,N 

KK  = A 

IF<INDEX(K)-KA)57,38,37 

37  CONTINUE 

38  ITEMPMNDEXU) 

INDEX(J)=INDEX(KK) 

INDEX(tCX)  = ITEMP 

700  CONTINUE 


IF(N)20P0. 1000, 2000 
2000  L *N - 1 


00  39  J * JA.NN 
LL  = 1 

00  *2  ( = 1,N 

C THIS  SECTION  IS  US  E 0 TO  SEE  IF  ONLY  A SINGULAR  TYPE  SOLUTION  IS  P 

IF  <CABS<A<1C  ,J  ) > - 0.0)  *3, *2. *3 
*2  CONTINUE 


I 7-  INOE  X <N> 

IF  <CA8S(A(IZ,N)  ) - 1.E-2)  *6,46,** 

**  CONTINUE 

WR1TE<3,*5> 

*5  FORMAT! 1 X , ONLY  SOLUTION  IS  ZERO  VECTOR*') 

60  TO  10 
*6  CONTINUE 

WR1 TE(3,*5) 

GO  TO  10 

c this  loop  is  equivalent  to  io.*,7  of  hiloebrand 


*3  DO  *0  I J =Lk.  ,l 
SUM1=0.0 
I I = N — I J 
I-INPEX(II) 
LLL  = 11*1 


DO  9300  K=LLL,N 
IP* INDE  X <X ) 

9300  SUM1*SUM1*A(I ,X) « A( IP  ,J ) 


b-i: 


171 

172 

173 

174 

175 

176 

177 

178 

179 

180 
181 
18? 

183 

184 

185 

186 

187 

188 

189 

190 

191 
19? 

193 

194 

195 

196 

197 

198 

199 

?00 


A ( 1 ,J>*A(l,J)-SUM1 
40  CONTINUE 

39  CONTINUE 


1000  CONTINUE 


C THIS  SECTION  SHIFTS  THE  SOLUTION  MATRIX  INTO  THE  FIRST  N BY  M 

C LOCATIONS  (ROUS  BY  COLUMNS) 


DO  400  1=1, N 
DO  400  J * JA.NN 
K=INOEX(I ) 

L = J-N 

400  A ( I ,L>  = A(K,J) 


COMMENT WRITE  ANSWER  INTO  X MATRIX 

DO  ?50  IQ5  = 1 , N 
DO  ?50  I 06  = 1 , M 

?50  X(I05,IQ6)  = A(IQ5«IQ6) 

10  CONTINUE 
RETURN 
END 
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SUBROUTINE  FOUR(DELCOS,AC.E,CI,CJ,NTERMS,NSAMPS> 

THIS  SUBROUTINE  CALCULATES  THE  FOURIER  C OE F F I C I ENT S ,C I AND  CJ.FOR  A 
FIELD  PATTERN, E, IN  EQUAL  STEPS  OF  COSINE  THETA. 

DIMENSION  AC(100),CI(51),CJ(50),E(100) 

DATA  PI/3.1415927/ 

CALCULATE  Cl  TER*S 
DO  20  1=1,NTERMS*1 
S UM  = 0 • 

N = I-1 

DO  10  J * 1 , N S AMPS 
AR6=N»PI*AC(J) 

SUM=SUM*E ( J ) * C 0 S ( AR6) 

CI(1)=DELC0S«SUM 

CI(T)=0.5*CIC1) 

CALCULATE  cj  terms 
DO  40  1 = 1 ,NTERMS 
SUM  = 0 • 

DO  30  J = 1 -NSAMPS 
AR6=1*P1*AC(J) 

SUM=SUM*E(J)*SIN( ARE) 

C J(  I)  = DELCOS*SUM 
RETURN 
END 
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SUBROUTINE  FOURO(D£LCOS .E.CI.NSAMPS) 

THIS  SUBROUTINE  CALCULATES  THE  FOURIER  COEFFICIENT  CKO)  FOR  A 
FIELD  PATTERN, E, IN  EQUAL  STEPS  OF  COSINE  THETA. 

DIMENSION  EC100) 

DATA  PI/3. 1415927/ 

CALCULATE  CKO)  TERM 
S UM  = Q , 

DO  10  J = 1 .NSAMPS 
) SUM  = SUM*  E ( J ) 

C I=0.5«0ELC0S*SUM 

RETURN 

END 


SUBROUTINE  MULL E R (C OE  ,C 0 E I ,N 1 . R 00 T R . RO OT I ) 
FORTRAN  MULLER  C2.2-001B 

MULLER 

DIMENSION  COE (51  ) .ROOTR ( 50  ),ROOT I (50 I.COEl <5 1 ) 
N2  = M ♦ 1 
N4-0 
I = NK1 

19  IF(C0E(I>)9,7.9 
7 I F C C OE 1(1))  9,1,9 
1 N4  = N4  ♦ 1 

R00TR(N4)=0, 

ROO  T I (N4 ) =0 • 

1 = 1-1 

I F ( N4  — N 1)19,37,19 
9 CONTINUE 

10  AXR-0.97 
AXI=0. 

N3=  1 

ALP1R=AXR 
ALP1 I=AXI 
M = 1 

60  TO  99 

11  BE  T 1 R = TE  MR 
BET  1 1MEM1 
AXR  =0.98 
ALP?R=AXR 
ALP  ? I *A  X 1 
M = ? 

60  TO  99 

12  BET2R=TEMR 
BET2I=TEMI 
AXR=0.99 
ALP3R=AXR 
ALP3I=AXI 
M = 3 

60  TO  99 

13  BET3R=TEMR 
BET  3I  = TEMI 

14  TE1 =ALP1 R-ALPJR 
TE2=ALP1 I -ALP3I 
TE5=ALP3R-ALP2R 
TE6=ALP3I -ALP2I 
TEM=TE5*TE5*TE6*TE6 
TE3=(TE1*TE5+TE2*TE6)/TEM 
TE4  = (TE2«TE5-T£1*TE6)  /TEN 
TE  7 = TE  3+  1 . 

TE9=TE3*TE3-TE4»TE4 

TE10=2.*TE3*TE4 

DEI 5 = TE7»BET3R-TE4*BE  T3I 
0E16  = TE7*BET3KTE4«BE  T3R 
TE11=TE3*BET2R-TE4*BET2I*BET1R-DE15 
TE12*TE3«BET2I*TE4»BET2R*8ET1I-DE16 
TE7=TE9-1 . 

TE1=TE9«BET2R-TE10*BET2I 
TE2  = TE9«BET2KTE10*PET2R 
TE13=T£1-8ET1R-TE7*BET3R*TE10»BFT31 
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TE14»TE2-BET1I-T£7«BE T5 1 -T E 1 0*8 E T 3* 

TE15XD£1J*TE3-DE16*TE4 

TE16x0E15»Te4*0E16*TE3 

TE1=TE13»TE13-TE14«TE14-4.*(TE11«TE15-TE12*TE16) 

TE2x2.*TE13«TE14-4.*(TE12*TE15*TE11*TE16> 

TENxS8RT(T£l»TE1*TE2*TE2) 

1MTEU113.113.112 
TE4=S0RT(.S*(TEN-TE1)  ) 

TE3=.5*T£2/TE4 
60  TO  111 

TE3XS9RT<  .5* (TEN*  TE1 ) ) 

IF(T£2)110,200,200 

TE3--TE3 

TE4x,5*TE2/TE3 

TE7=TE13*TE3 

TE8=TE14*TE4 

TE9=TE13-TE3 

TE10xT£l4-TE4 

TE1=2.«T£15 

TE2=2.*TE16 

IE ( TE7*TE7*TE8*TE8-TE  9«TE9-TE10«TE10  )204  ,204  ,205 

TE7=TE9 

TE8=TE10 

TENXTE7»TE7*TE8»TE8 

TE3=(TE1«TE7*TE2*TE8)/TEN 

TE4=(TE2*TE7-TE1*TE8)/TEN 

AXRxALP3R*TE3»TE5-TE4»TE6 

AXI xALP3I*TE3«TE6*TE4*TE5 

ALP4R-AXR 

ALP4I-AXI 

N=4 

60  TO  99 
N6X  1 

IE  ( ABS(  he  LU  + ABS  < BELL  )-1  ,E -20)  18, 18, 16 

TE7XABS( ALP3R-AXR )*ABS( ALP3I-AX!) 

IF(TE7/(ABS(AXR)*ABS(AXI>)-1.E-t>18,18,17 

N3XN3*1 

ALP  1RXALP2R 

ALP  1 1*ALP2I 

ALP2RxALP3R 

ALP21xALP31 

ALP3R=ALP4R 

ALP3l*ALP4I 

BET1RXBET2R 

BET1I=BET2l 

8ET2R*BE T3R 

BET2IxB£T3l 

BET3R=TENR 

8ET31XTENI 

I F ( N3  .LT.  100)  GO  TO  14 

I FL  A6*1 

WRI TE (3,400> 

EORNATMX,  ‘--ERROR  — BAD  ROOTS  IN  SUBROUTINE  PULLER*') 
N4*  N4 ♦ 1 

R00TR(N4)-ALP4R 
ROOTI (N4)xALP4I 

I F ( I F LA6  .6  9.  1)  URITE(3,4Q1)  N4  ,R00 TR <N 4 ) ,R 00 T I (N 4 ) 


11* 

m 

117 
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119 
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Ml 

Mi 

126 

127 

128 
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w 

132 

133 

Mi 

136 


*01  F0I)NAT(1X,I3,2E15.S) 

IELAS=0 
N3*  0 

*1  If <N*-N1 110,37,37 
37  RETURN 
99  TENR=C0E(1> 

ten  1=0.0 

00  100  1=1, HI 

TE1 =TENR»AXR-T£NI »AXI 
TENI=TENI*AXR*TENR«AXI  ♦COEKI  + 1) 
100  TEN  R = TE1«C0E(I«1) 

MELL=T£NR 

BELL=TENI 

*2  IF(N*)102,103,102 

102  00  101  1=1, N* 

TENlsAXR-ROOTRU  ) 

TEN?=AXI -ROOTI (1 ) 
TE1=TEN1«TEN1«TEN2*TEN2 
TE2=<T£NR*TEN1*TENI«TEN2)/TE1 
TENI*(TENI«TEN1-TENR«TEN2)/TE1 


m 


TEN  R = TE  2 

0(11,  12, 13, 15), N 


«0  TOC 
END 
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SUBROUTINE  PRON Y ( DEL  X ,Y V AL $ . A , ALPHA , NT ER NS ,N PO IN T ) 
OINENSION  RINRC51  >»  RINM51),  ROOTR(S0>,  R00TK50) 
CONPLEX  YVALS,  f,  fBARF,  B,  FBARB , M,  F1BARY.EIG, 

1 SOLN,  A,  ALPHA, F1BARF 

C OH  PON  /BLOCKl/  F < 1 00  ,5 0 ) , FB AR F < 50 , 5 0 ) ,B (1 00 ) . f B ARB ( 5 0 ) 
EOU I VALENCE  <f. FI), <F BA RF.F1BARF), (FBARB, FIBARY) 
DIMENSION  F 1 M 00 , 50  > , F 1 B AR F< 5 0 , 50 ) . F1BART(50) 
DIMENSION  SOLN (50  ),  A (50),  AlPHA(50) 

DIMENSION  T VAL  S C 1 00 ) 

1DIFF  = NPOINT  - NT  E R MS 


430 

200 


If ( IDIFF 
WR1TE(3, 43 


128)  GO  TO  200 


• LE.  1 

■ M • C \ J t *30  ) 
f OR  MAT  (1  X,’* --ERR  OR—  NPTS-NPOLES 

CONTINUE 


.GT.  128*') 


COMMENT- 


--NON  TO  GENERATE  F MATRIX 
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DO  800  IA*1,  NTERHS 

00  600  16=1,  IDIFF 

IPASS=NTERHS 

800  F(IB.IA)  * YVALS ( IA« I B-1 ) 

COMMENT NOW  GENERATE  B VECTOR 

DO  801  I E = 1 , IDIFF 

801  B ( I E ) = - YVALS(IE*IPASS) 

COMMENT COMPUTE  FBARF  AND  FBARB 

DO  802  1*1, NTERMS 
DO  802  J=1,NTERMS 
FBARF (I 

DO  802  K=1, IDIFF 

802  FBARF (I ,J)=FBARF (I , J )*F (K  ,1 )*F ( K , J ) 

DO  803  1=1, NTERHS 
FB A RB ( I )=0, 

DO  803  K*1, IDIFF 

803  FBARB(I)=FBARB(1)*F(X,I)*B(K) 

COMMENT USE  CROUT  TO  SOLVE  EON  f B A R F *S = FB AR B FOR  S 

CALL  CROUMFBARF  , SO  L N , FBARB,  N T E RMS  , 1) 

COMMENT CIN  NON  CONTAINS  SOLUTIONS  FOR  S 

COMMENT NOW  TO  FIND  ROOTS  OF  POLYNOMIAL  WITH  S VECTOR  AS  COEFFS 

RINR(1>=1 . 

RINIC1)=0, 

DO  604  I M*  1 .NTERMS 
IZ*NTE RMS*1-IM 
I Y= I M* 1 

RINR(IY)  = REAL(S0LN(1Z)) 

804  RINI(IY)  * A IMAG ( SOLN  (I Z ) ) 

CALL  MULLER  (RINR,  RINI,  NTERMS,  ROOTR,  ROOTI) 


